{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Storage of Network Data and Topology\n",
    "\n",
    "As discussed in the first tutorial, OpenPNM uses Numpy ``ndarrays`` to store all data, and these arrays are all stored in ``dict``s, because all OpenPNM objects are subclassed python dictionaries. In this tutorial will take a deeper dive into the details of the OpenPNM data storage scheme. The sections to be covered are:\n",
    "\n",
    "- The Spreadsheet Analogy \n",
    "- Maintaining Data Integrity\n",
    "- Representing Topology"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Spreadsheet Analogy\n",
    "\n",
    "The best analogy for explaining data storage in OpenPNM is the humble spreadsheet.  According to this analogy you can imagine all data is stored in a table; two tables actually, one for pore data and one for throat data.  Each pore (or throat) corresponds to a row and each property corresponds to a column.  Consider the following network with 4 pores, 3 throats:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "import openpnm as op\n",
    "import numpy as np\n",
    "op.visualization.set_mpl_style()\n",
    "np.random.seed(0)\n",
    "pn = op.network.Demo(shape=[4, 1, 1])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's use [pandas](https://pandas.pydata.org/) to express the geometric properties as a 'spreadsheet' and pull out all pore properties that are as wide as a single column."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "pore_data_sheet = pd.DataFrame({k: pn[k] for k in pn.props(element='pore') if pn[k].ndim == 1})"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can now view this 'spreadsheet':"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "application/vnd.microsoft.datawrangler.viewer.v0+json": {
       "columns": [
        {
         "name": "index",
         "rawType": "int64",
         "type": "integer"
        },
        {
         "name": "pore.coordination_number",
         "rawType": "float64",
         "type": "float"
        },
        {
         "name": "pore.diameter",
         "rawType": "float64",
         "type": "float"
        },
        {
         "name": "pore.max_size",
         "rawType": "float64",
         "type": "float"
        },
        {
         "name": "pore.seed",
         "rawType": "float64",
         "type": "float"
        },
        {
         "name": "pore.volume",
         "rawType": "float64",
         "type": "float"
        }
       ],
       "ref": "82c7bfa6-1c6d-46db-bdf0-8c3cbd8b7619",
       "rows": [
        [
         "0",
         "1.0",
         "0.47440675196366233",
         "1.0",
         "0.47440675196366233",
         "0.05590507143096387"
        ],
        [
         "1",
         "2.0",
         "0.5575946831862097",
         "1.0",
         "0.5575946831862097",
         "0.09077254130465309"
        ],
        [
         "2",
         "2.0",
         "0.501381688035822",
         "1.0",
         "0.501381688035822",
         "0.06599393532741361"
        ],
        [
         "3",
         "1.0",
         "0.4724415914984484",
         "1.0",
         "0.4724415914984484",
         "0.05521320959471004"
        ]
       ],
       "shape": {
        "columns": 5,
        "rows": 4
       }
      },
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>pore.coordination_number</th>\n",
       "      <th>pore.diameter</th>\n",
       "      <th>pore.max_size</th>\n",
       "      <th>pore.seed</th>\n",
       "      <th>pore.volume</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>1.0</td>\n",
       "      <td>0.474407</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.474407</td>\n",
       "      <td>0.055905</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2.0</td>\n",
       "      <td>0.557595</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.557595</td>\n",
       "      <td>0.090773</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2.0</td>\n",
       "      <td>0.501382</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.501382</td>\n",
       "      <td>0.065994</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>1.0</td>\n",
       "      <td>0.472442</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0.472442</td>\n",
       "      <td>0.055213</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   pore.coordination_number  pore.diameter  pore.max_size  pore.seed  \\\n",
       "0                       1.0       0.474407            1.0   0.474407   \n",
       "1                       2.0       0.557595            1.0   0.557595   \n",
       "2                       2.0       0.501382            1.0   0.501382   \n",
       "3                       1.0       0.472442            1.0   0.472442   \n",
       "\n",
       "   pore.volume  \n",
       "0     0.055905  \n",
       "1     0.090773  \n",
       "2     0.065994  \n",
       "3     0.055213  "
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pore_data_sheet"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The column names are the properties such as **pore.volume**, and the rows correspond to the pore index, so pore ``0`` has a volume of 0.055905.  "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "One could also extract an entire column using:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-06-24T11:25:06.060457Z",
     "iopub.status.busy": "2021-06-24T11:25:06.059315Z",
     "iopub.status.idle": "2021-06-24T11:25:06.063503Z",
     "shell.execute_reply": "2021-06-24T11:25:06.064416Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0    0.055905\n",
      "1    0.090773\n",
      "2    0.065994\n",
      "3    0.055213\n",
      "Name: pore.volume, dtype: float64\n"
     ]
    }
   ],
   "source": [
    "column = pore_data_sheet['pore.volume']\n",
    "print(column)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Or access individual elements:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-06-24T11:25:06.071219Z",
     "iopub.status.busy": "2021-06-24T11:25:06.070328Z",
     "iopub.status.idle": "2021-06-24T11:25:06.075321Z",
     "shell.execute_reply": "2021-06-24T11:25:06.074525Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.05590507143096387\n"
     ]
    }
   ],
   "source": [
    "print(pore_data_sheet['pore.volume'][0])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "```{warning}\n",
    "The spreadsheet analogy is very apt, but it breaks down slightly when one considers data that is multiple columns wide such as **'pore.coords'**.  The is stored as an Np-by-3 array in the network dictionary, but spreadsheets usually require that all data is only a single column wide. \n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Rules to Maintain Data Integrity\n",
    "\n",
    "Several rules have been implemented to control the integrity of the data. This is implemented by subclassing the `__setitem__` method of the `dict`. Each time data is assigned to a dictionary key, OpenPNM first checks the data type and dictionary key to ensure the rules are followed.  Each of these checks are explored below.\n",
    "\n",
    "#### All Values are Converted to Numpy Arrays\n",
    "Only Numpy arrays can be stored in an OpenPNM object, and any data that is written to an OpenPNM object (i.e. dictionary) will be converted to a Numpy array.  This is done to ensure that all mathematical operations throughout the code can be consistently done using vectorization. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-06-24T11:25:06.176510Z",
     "iopub.status.busy": "2021-06-24T11:25:06.175396Z",
     "iopub.status.idle": "2021-06-24T11:25:06.180431Z",
     "shell.execute_reply": "2021-06-24T11:25:06.181503Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "<class 'numpy.ndarray'>\n"
     ]
    }
   ],
   "source": [
    "pn['throat.list'] = [1, 2, 3]\n",
    "print(type(pn['throat.list']))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This illustrates that the basic python list-type has been converted to a Numpy array when stored in the dictionary"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Dictionary Keys Must Start With 'pore' or 'throat'\n",
    "All array names must begin with either *'pore.'* or *'throat.'* which serves to identify the type of information they contain."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-06-24T11:25:06.204358Z",
     "iopub.status.busy": "2021-06-24T11:25:06.203478Z",
     "iopub.status.idle": "2021-06-24T11:25:06.208127Z",
     "shell.execute_reply": "2021-06-24T11:25:06.208790Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "This will throw an exception since the dict name cannot start with foo\n"
     ]
    }
   ],
   "source": [
    "try:\n",
    "    pn['foo.bar'] = 0\n",
    "except:  # noqa E722\n",
    "    print('This will throw an exception since the dict name cannot start with foo')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The pore or throat prefix is enforced so that OpenPNM knows how long an array should be. Attempting to write an array of the wrong length will result in an error:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "This will throw an exception since there are 4 pores, hence all pore arrays should be 4 elements long\n"
     ]
    }
   ],
   "source": [
    "try:\n",
    "    pn['pore.test'] = [0, 0, 0]\n",
    "except:  # noqa E722\n",
    "    print('This will throw an exception since there are 4 pores, hence all pore arrays should be 4 elements long')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Any Scalars are Expanded to a Full-Length Vector\n",
    "For the sake of consistency only arrays of length *Np* or *Nt* are allowed in the dictionary. Assigning a scalar value to a dictionary results in the creation of a full length vector, either *Np* or *Nt* long, depending on the name of the array. This effectively applies the scalar value to all locations in the network."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-06-24T11:25:06.191724Z",
     "iopub.status.busy": "2021-06-24T11:25:06.190288Z",
     "iopub.status.idle": "2021-06-24T11:25:06.195032Z",
     "shell.execute_reply": "2021-06-24T11:25:06.196182Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[0 0 0 0]\n"
     ]
    }
   ],
   "source": [
    "pn['pore.test'] = 0\n",
    "print(pn['pore.test'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Note how the scalar value has been cast to an array of 4 elements long, one for each pore in the network. "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Nested Dictionary Names are Allowed\n",
    "It's possible to create nested properties by assigning a dictionary containing numpy arrays"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-06-24T11:25:06.218609Z",
     "iopub.status.busy": "2021-06-24T11:25:06.217483Z",
     "iopub.status.idle": "2021-06-24T11:25:06.222264Z",
     "shell.execute_reply": "2021-06-24T11:25:06.223411Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{'species_A': array([0, 0, 0, 0]), 'species_B': array([1, 1, 1, 1])}\n"
     ]
    }
   ],
   "source": [
    "pn['pore.concentration'] = {'species_A': 0, 'species_B': 1}\n",
    "print(pn['pore.concentration'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The above rule about expanding the scalar values to a numpy array have been applied.  \n",
    "\n",
    "Requesting the top level of dictionary key returns both concentrations, but they can be accessed directly too:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-06-24T11:25:06.231472Z",
     "iopub.status.busy": "2021-06-24T11:25:06.230502Z",
     "iopub.status.idle": "2021-06-24T11:25:06.235503Z",
     "shell.execute_reply": "2021-06-24T11:25:06.236221Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[0 0 0 0]\n"
     ]
    }
   ],
   "source": [
    "print(pn['pore.concentration.species_A'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "You can also retrieve a dictionary by requesting only the element ('pore' or 'throat') and the propname ('concentration'):"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-06-24T11:25:06.244270Z",
     "iopub.status.busy": "2021-06-24T11:25:06.243340Z",
     "iopub.status.idle": "2021-06-24T11:25:06.248217Z",
     "shell.execute_reply": "2021-06-24T11:25:06.249026Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'species_A': array([0, 0, 0, 0]), 'species_B': array([1, 1, 1, 1])}"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pn['pore.concentration']"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The above returns a 'subdictionary' which can be indexed into:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([0, 0, 0, 0])"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pn['pore.concentration']['species_A']"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Boolean Arrays are Treated as Labels, Numerical Arrays are Treated as Properties\n",
    "Any Boolean data will be treated as a *label* while all other numerical data is treated as a *property*. OpenPNM uses *labels* extensively as will be illustrated in a later tutorial. For this tutorial, we'll just introduce the idea.  We can create a new label by assigning `False` to a new key. This will expand the single `False` value into every pore, which means that NO pores have that label (yet):"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-06-24T11:25:06.257992Z",
     "iopub.status.busy": "2021-06-24T11:25:06.256769Z",
     "iopub.status.idle": "2021-06-24T11:25:06.261857Z",
     "shell.execute_reply": "2021-06-24T11:25:06.262995Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "['pore.label', 'pore.left', 'pore.right', 'pore.surface', 'pore.xmax', 'pore.xmin']\n"
     ]
    }
   ],
   "source": [
    "pn['pore.label'] = False\n",
    "print(pn.labels(element='pore'))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "You can see that ``'pore.label'`` shows up in this list automatically since it is of Boolean type. However, if you ask the network which pores have that label, the answer is none:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([], dtype=int64)"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pn.pores('label')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "If we then set some locations to `True`, this will change:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[0 1 2]\n"
     ]
    }
   ],
   "source": [
    "pn['pore.label'][[0, 1, 2]] = True\n",
    "print(pn.pores('label'))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "OpenPNM objects include a method for creating a label and applying it to locations all in one step using the `set_label` method."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[2 3]\n"
     ]
    }
   ],
   "source": [
    "pn.set_label('another', pores=[2, 3])\n",
    "print(pn.pores('another'))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Dictionary Keys with a Leading Underscore are Hidden\n",
    "Following the Python convention, if a piece of data is not really meant to be seen or used by the user, it can be pre-pended with an underscore and it will not appear in any output."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-06-24T11:25:06.271759Z",
     "iopub.status.busy": "2021-06-24T11:25:06.270788Z",
     "iopub.status.idle": "2021-06-24T11:25:06.275494Z",
     "shell.execute_reply": "2021-06-24T11:25:06.276273Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――\n",
      "1     : pore.concentration.species_A\n",
      "2     : pore.concentration.species_B\n",
      "3     : pore.coordination_number\n",
      "4     : pore.coords\n",
      "5     : pore.diameter\n",
      "6     : pore.max_size\n",
      "7     : pore.seed\n",
      "8     : pore.test\n",
      "9     : pore.volume\n",
      "10    : throat.conns\n",
      "11    : throat.cross_sectional_area\n",
      "12    : throat.diameter\n",
      "13    : throat.diffusive_size_factors\n",
      "14    : throat.hydraulic_size_factors\n",
      "15    : throat.length\n",
      "16    : throat.lens_volume\n",
      "17    : throat.list\n",
      "18    : throat.max_size\n",
      "19    : throat.spacing\n",
      "20    : throat.total_volume\n",
      "21    : throat.volume\n",
      "――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――\n"
     ]
    }
   ],
   "source": [
    "pn['pore._hidden'] = 1\n",
    "print(pn.props())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The ``'pore._hidden'`` key does not show up in this list, though it is present:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([1, 1, 1, 1])"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pn['pore._hidden']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "True\n"
     ]
    }
   ],
   "source": [
    "print('pore._hidden' in pn.keys())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## The `param` prefix\n",
    "\n",
    "Expanding all scalar values to arrays that are either Np or Nt long is done as a convenience or shorthand to mean \"place this value of X in all pores\". It means that you can later ask \"what is the value of X in pore Y\" and get an answer. This implies that X might take on different values in each pore/throat. However, there are properties which we know are constant, such as the molecular weight or critical temperature of a pure species. In V3 we have introduced the `params` attribute to all object, such that `pn.params` is a dictionary containing such scalar values. For convenience we have also added the ability to read and write to the `params` dictionary via the main get/set method. This is illustrated below:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――\n",
      "Parameters                          Value\n",
      "――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――\n",
      "dimensionality                      (3,)\n",
      "lattice_connectivity                6\n",
      "――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――\n"
     ]
    }
   ],
   "source": [
    "pn.params['lattice_connectivity'] = 6\n",
    "print(pn.params)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "```{tip} \n",
    "  **Use the \"param\" prefix:** It is possible use the `param` prefix when reading and writing to the object's dictionary. This intercepts the given key and dispatches the request to the `params` attribute. This feature is particularly useful when passing arguments to pore-scale models (to be discussed in a later tutorial) since you can pass either `arg='param.value'` or `arg='pore.value'` and OpenPNM will be able to fetch the requested data by doing `value = pn[arg]` and both options will work.\n",
    "```"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――\n",
      "Parameters                          Value\n",
      "――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――\n",
      "dimensionality                      (3,)\n",
      "lattice_connectivity                6\n",
      "test                                2\n",
      "――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――――\n"
     ]
    }
   ],
   "source": [
    "pn['param.test'] = 2\n",
    "print(pn.params)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Representing Topology\n",
    "\n",
    "Topology was introduced in the last tutorial, but not discussed in depth. Actually, the *next* tutorial will explain topology in full detail, but here we can reveal a little more about it as it pertains to data storage."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Consider the following simple random network:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-06-24T11:25:06.285814Z",
     "iopub.status.busy": "2021-06-24T11:25:06.284955Z",
     "iopub.status.idle": "2021-06-24T11:25:06.303321Z",
     "shell.execute_reply": "2021-06-24T11:25:06.304283Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeIAAAUxCAYAAABUDIUxAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjcsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvTLEjVAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAj1pJREFUeJzt3Qm4jeX6x/GfmSgiQqgcGlHRrEHznAblVFIqhJA6p/5NpzqdTtMpIWSIpJkmmktzGgylUkIjEZEpU4r/de/H6l3Lntbee631rPW+3891uXqWYe+nd7PvdT/DfZfbtGnTJgEAAC/K+/m0AADAEIgBAPCIQAwAgEcEYgAAPCIQAwDgEYEYAACPCMQAAHhEIAYAwCMCMQAAHhGIAQDwiEAMAIBHBGIAADwiEAMA4BGBGAAAjwjEAAB4RCAGAMAjAjEAAB4RiAEA8IhADACARwRiAAA8IhADAOARgRgAAI8IxAAAeEQgBgDAIwIxAAAeEYgBAPCIQAwAgEcEYgAAPCIQAwDgEYEYAACPCMQAAHhEIAYAwCMCMQAAHhGIAQDwiEAMAIBHBGIAADwiEAMA4BGBGAAAjwjEAAB4RCAGAMAjAjEAAB4RiAEA8IhADACARwRiAAA8IhADAOARgRgAAI8IxAAAeEQgBgDAIwIxAAAeEYgBAPCIQAwAgEcEYgAAPCIQAwDgEYEYAACPCMQAAHhEIAYAwCMCMQAAHhGIAQDwiEAMAIBHBGIAADwiEAMA4BGBGAAAjwjEAAB4RCAGAMAjAjEAAB4RiAEA8IhADACARwRiAAA8IhADAOARgRgAAI8IxAAAeEQgBgDAIwIxAAAeEYgBAPCIQAwAgEcEYgAAPCIQAwDgEYEYAACPCMQAAHhEIAYAwCMCMQAAHhGIAQDwiEAMAIBHBGIAADwiEAMA4BGBGAAAjwjEAAB4RCAGAMAjAjEAAB4RiAEA8IhADACARwRiAAA8IhADAOARgRgAAI8IxAAAeEQgBgDAIwIxAAAeEYgBAPCIQAwAgEcEYgAAPCIQAwDgEYEYAACPCMQAAHhEIAYAwCMCMQAAHhGIAQDwiEAMAIBHBGIAADwiEAMA4BGBGAAAjwjEAAB4RCAGAMAjAjEAAB4RiAEA8IhADACARwRiAAA8IhADAOARgRgAAI8IxAAAeEQgBgDAIwIxAAAeEYgBAPCIQAwAgEcEYgAAPCIQAwDgEYEYAACPCMQAAHhEIAYAwCMCMQAAHhGIAQDwiEAMAIBHBGIAADwiEAMA4BGBGAAAjwjEAAB4RCAGAMAjAjEAAB4RiAEA8IhADACARwRiAAA8IhADAOARgRgAAI8IxAAAeEQgBgDAIwIxAAAeEYgBAPCIQAwAgEcEYgAAPCIQAwDgEYEYAACPKvr85EAqbNokLVggTZ8uLV4s/f67VLmyVK+e1Lq11LChVK4czxpAdiIQIydt3Ci9/ro0YoT07rvSokWF/97tt5cOPVTq2lU6+mipPOtAALJIuU2bLJ8AcsOff0rDh0v33CPNnVvyP9+smXTFFVK3blKFCumYIQCUDIEYOWPOHKlLF+n99xN/vlYtad99pTZtpL/9zS1L2/L0N99I06ZJU6dKy5cn/pm2baXRo6XmzTP6vwAA+RCIkRPGjnVZ7Lp1wc8ddZTUs6d06qlSxSI2Wf74Q5owQRoyRJo0Kfj5qlVddn3++emdOwAUhUCMrGcBtFev4LVlvaNGSYcdVvKP9c47Lqv+9tvg5wYPdgEdAHzg2AqyPhOOD8J24GrGjNIFYWN/7rPP3MeJsY9vnwcAfCAjRlbvCbdqFSxHX3ONdOutqbmKZEcUr7tOuu0297paNRfg2TMGkGkEYmTt6ejDDw8OZlkGO2xYau8DWzDu3t1dgYod4Hr7bU5TA8gsAjGy0tChwb6t7Qlbtlq9euo/z+rVLuuO7RnbfnSPHqn/PABQGPaIkZXFOuyecIwdzCoqCE+ZMkWXXXaZ9txzT1WvXl1NmjTR2WefrdmzZxf7uezj2seP6d/ffX4AyBQyYmSdV1+VjjsuuKJkFbSK0qFDB73//vs666yz1KpVK/3888+677779Ntvv+nDDz9UixYtiv2c9nneeCP4/Mcck4r/EwAoHoEYWeess6Tx4934qaekM84o+vdPnjxZ++67rypbJY/N5syZo5YtW+YF6YcffrjYz2mfp0MHN7b/jhtXtv8HAEgWgRhZxQ5QNWjgakdbxaxffim6WEdR2lipLVl1rWnF/l4r+lG3rqvAVb++tHBh6T4nAJQUe8TIKtZFKdbAwcpWljYIWwn1RYsWabvttkvq99vn2Ry39fPPbh4AkAkEYmQVa2UYEwuMpfHII4/op59+UseOHZP+M/GfL4kkGgBSgkCMrGL9hGPs2lJpzJo1S7169dJBBx2kCy64IOk/F//5bEkcADKBQIysYl2TYuLOXiXNTkyfdNJJqlmzpsaPH68KJeh1WKVKMF6/vuSfGwBKo5Q7cEB6xAff+KCcjBUrVuiEE07Q8uXL9e6776phw4Yl+vPxwTc+KANAOhGIkVXq1QvG1k84WevWrdMpp5ySV8Tj9ddf1x577FHizx3/+ewENQBkAoEYWaV165IfmPrzzz/zDmV98MEHeu655/L2hksj/vOV5aAYAJQE94iR8/eIL7/8cg0YMCAvI7bSllvq1KlTsZ+Xe8QAfCEQI+cra7Vr105vW9ukIu4UF4fKWgB84dQ0so61PIyxbkjFeeutt/KCbWE/khH/ebp1K82sAaB0yIiRdaz70a67SnPnuteW7B52WPo+n338du3cuHlzu4cslectKoAM4dsNso4FwSuuCF5fdJHrG5wO9nHt48f060cQBpBZBGJkJVsebts2uFZkATLJVeak2ce7/HLp22/d60MOkbp3T+3nAIDiEIiRlawg1ujRUtWq7vWIEdJ116UuGNvHufZaaeRI97paNWnUKLJhAJlHIEbWsv3a4cOD17fd5jLWsi5T25+3jPv224OfGzbMfT4AyDQOayHr2YnmXr2C102buuz18MNLdzDL9oRjy9Fm8GCpZ8/UzBUASopAjJwwdqxlw5u0dm25v37uyCNdAG3fvuiiHxs2SBMmuID+xhvBz9tytGXC55+f5skDQBEIxMgZc24bry7X1tf7OiTh560Cl5WktB/WytAaNlgDBzvkZWUr7cfy5Ykfq63e0+ixldS80wGZ/Z8AgC0QiJEb1qyRdtlFf/60UMPVTffscI/m/lStxB+mWd3luuKXa9RNw1Vhl2bSjBnBiTAA8IDDWsgNgwZJP/2kCtqoHqf8pK9/rKZXX5U6dJDq1y/6j9qv2++z3//1gm3U46AZeR9Hs2e7E2AA4BEZMbLf0qVuzXnFCne/6LPPpD33TPgtCxa4JWhrEmHL0rY8ba0Mbbk6X1viL76Q9tnHdXqoVMllxbvvntH/JQCIoQ0isp9lrRaEzYUX5gvCxoJtvoBbmBYtpH/+031cO8lld5nsODV1LQF4QEaM7PbDD3l7w/r9d7eXO2eO1KhR2T/u2rVSy5buRJexC8vx3SYAIEPYI0Z2+9e/XBA2ffumJgjH7i7df3/w+qqrpJ9/Ts3HBoASICNG9rK9W9vLtXqU227rslf7byp17uwuKZu//1167LHUfnwAKAYZMbLXNdcExaWt0HSqg7C5+26pTh03fvxx6aWXUv85AKAIZMTITm++6UpnmSZNpK+/Tt993wcflLp0ceMdd5RmzpSqV0/P5wKALZARI/tYFmx7tjG33JLeohsXXCAdcURwOOymm9L3uQBgC2TEyD5PPil17OjGrVpJ06e7vojpZKex7RS1XUK2zzVlitufBoA0IyNGdrET0tYoOMZ6FaY7CBvrgXj99W7855/ubrH9FwDSjECM7DJiRHC315aLjz8+c5/blsP32MONp06V7rsvc58bQGSxNI3ssWqVK2VpdSrNxx9L++2X2Tm8/750yObuTjVqSF9+KTVunNk5AIgUMmJkD7tKFAvCZ5+d+SBs2ra1xsdu/NtvUq9ewRUqAEgDMmJkB6tq1ayZtHq1VLGi9NVX7rUP1rzYmkDEKm2NHy+deaafuQAIPTJiZAe7omRB2FhG6isIm1q1pAEDgte9ewdNJwAgxciI4Z9dHbJDUtaW0PZl586Vtt/e75xsOfqUU6QXXnCve/aUBg/2OycAoURGDP+sfKUFYfOPf/gPwqZcORd4t9rKvR46VPrgA9+zAhBCZMTwy05GH3CAG9er57LhrbfOnq/KPfdIV14Z9DG24iKVKvmeFYAQISNG9pSyvPHG7ArCpk8fqXVrN/7iC+l///M9IwAhQ0YMf158UTrpJDe2w1l2Zzcbs81p06T995c2bnQ1ry0g231nAEgBMmL4YeUjr746eP3f/2ZnEDZt2kh9+7rxunXSpZdytxhAypARw48xY6QLL3RjyzY//NAdkMpWVtxjzz2lH390rx96SDr/fN+zAhACBGJknmWVu+wizZsX9B5u1y77vxLPP++uNJnttnNFR+y/AFAGLE0j86yZQiwIn3hibgRhc/LJ0llnufGSJdI//+l7RgBCgIwYmbVsmTvoZP+1pehPP3U9h3PFwoWu/GWs0takSdKRR/qeFYAcRkaMzLL+whaETefOuRWETYMG7v8hxg5u2VI7AJQSGTEyx5ajmzeX1q+XqlSRZs+WmjTJva+AXWM69FBp8mT3+vrrXa1sACgFMmJkjhXssCAca6SQi0HYlC8vDR8eXLe64w53BxoASoGMGJlhRTD22stlk9bd6JtvpNq1c/vpWyZ8661BH+N33nFBGgBKgO8ayIxrrnFBODbO9SAca1YRa9f4/vvSiBG+ZwQgB5ERI/0sUzz8cDdu1MjtDVerFo4nb6emjz7ajWvWdHeL7UAXACSJjBiZbezw73+HJwibo46SLrjAje1K0+WX+54RgBxDRoz0euopqUMHN7YSkTNmSBUqhOupW3GP3XaTli4NKnDFmlkAQDHIiJE+Gza4/eAYu38btiBsrMyl9S2O6dnT1aYGgCQQiJE+DzwgzZnjxocdFu4s0RpA2DK1scYQdlULAJLA0jTSwzJCO1G8aJF7/cEH0oEHhvtpz50rtWzpKm3ZNaYpU6TWrX3PCkCWIyNGevTvHwThM88MfxA29sbjhhvc2K5qdesm/fGH71kByHJkxEi9xYtdYwfLim1PeOZMadddo/Gkf//dZcH2/2xs77hfP9+zApDFyIiRev/5T3BYqWvX6ARhU7myK38ZYxmy7RkDQCHIiJFaVrrS2gTaiemttnKv69eP3lO2k9NDhwZ9jCdMcG0fAWALZMRIff1lC8LmyiujGYTNbbcFFbbsXrHdpwaAApARI3WmTpX228+N69Z1p4i32Sa6T3j8eOmss9zY3pBY+UtreAEAcciIkbpSlldfnbg3GuUgHDstbsvS5uefE4ubAMBmZMRIjVdekY4/3o2bNnXZnx1cijo7qLXHHtLq1UGXpoMP9j0rAFmEjBhlZ3dm47Nh69FLEHaaNHGnyGPsbrFdcQKAzQjEKLtHH3XNHEybNtLZZ/NU4/Xu7Z6LsfvFd93F8wHwF5amUTZWztE6D/3wg3v9+utBzWUEpk93B9ls9aBKFenzz6XmzXlCAMiIUUZ2VzYWhI87jiBcGKu2FauwtX69dOml7oAbgMgjI0bpLV/uSln++qsrVmFZ395780QLYwe2rCdz7I3LmDFS5848LyDi2CNG6d15pwvC5rzzCMLFqV5dGjIkeH3FFdKSJfwNBCKOjBil89NPbo9z7Vp3Qvrrr6WdduJpJqNjR+nJJ93YMmLLjAFEFhkxSuemm1wQNr16EYRLYsAAqWZNN37oIWnSJP4WAhFGRoyS+/JLqWVLdwLYqmd9+61Upw5PsiSsQ1P37kEf488+k6pV4xkCEURGjJKzUo0WhM3//R9BuDQuuURq29aNrSZ3fNEPAJFCRoySee896dBD3bhhQ2nOHNfuEKVbWbBT5tatqmJF6ZNPpBYteJJAxJARI3l27/Wqq4LXN99MEC4Lq0FtKwrmjz/cUnVspQFAZJARI3nPPiudfrob776729e0TA5lq0y2117S7NlBgRQr9gEgMgjESI5lbHZAa9asICi3b8/TS4U335SOPNKN7fCbda6yZX8AkcDSNJIzenQQhO2Q0amn8uRS5YgjpAsvdOOVK6W+fXm2QISQESO50oxWvGPhQveanrqpt3Spa54Rq7Q1caJ08sn87QQigIwYyRWgiAXh006jsX062D3s/v2D11Yk5bff+NsJRAAZMYpmGZo1drAl0/LlXT9dy9yQnlPp1sHqtdfca+vWdM89PGkg5MiIUbRbb3VB2Fx8MUE4nayDlZ2arlo1WImYNo2/oUDIkRGjcN99J+26qys4YeUXrQIUp3nT7/bbXfUys88+0scfc00MCDEyYhTu+utdEI4tkxKEM+PKK91VMWPVtgYOzNAnBuADGTEKZgGgdevgINE33wQdg5B+H37oDsXZvrGVELW9edpMAqFERoyCXX11YmZMEM6sAw+UevRw4zVr3ClqC8oAQoeMGPnZqd1jj3Vjy8KskEeVKjypTFuxwtWjXrDAvX7iCenss/k6ACFDRoxE1nQgPhu29nwEYT9sFWLQoOC1VdxavtzTZACkC4EYiSzrsv3h2Indc87hCflkTTZi5UR//jno1gQgNFiaRmD9etdVya4tmVdflY45hifk27x5bok6Vmnr3XelQw7xPSsAKUJGjMCwYUEQPvpognC2aNzYFVaJsb7Fv//uc0YAUoiMGI5Vz7JSlrGmA1bRKXZ9Cf79+ad00EHSlCnu9S23uNPsAHIeGTGcu+4KgrDtCxOEs0uFCtLw4e6/sUN0s2f7nhWAFCAjhuus1KyZu69aqZK7rtS0KU8mG/3zn9L//hf0MZ40ydWoBpCzyIgh3XSTC8LGikgQhLP7axWrsPXmm9KYMb5nBKCMyIijzrLfFi3cHuTWW7tSlnXr+p4VivLyy9IJJ7hx7drua8jXDMhZZMRRd+21Lgibq67iG3ouOP546e9/d+Nff3VNIgDkLDLiKPvgA9dYwNSv79ocVq/ue1ZIxqJFrjd0rNIWd76BnEVGHGU1arjmArG9R4Jw7th+e3fSPcb29teu9TkjAKVERhxl1mvYTklbkwc7gVuxou8ZoaR1wdu1c5W2zDXXSP/9L88QyDEEYrhv6OVZHMlJX30l7bWXe1Nlb6SmT5datvQ9KwAlwHdfEIRzmdUGt0zY/PGH1K2be2MFIGeQEYfYxo0bVa5cubwfCLF166S995a+/tq9HjxY6tnT96wAJImMOITGjRun7777TuXLl88LwhaQEWJVq7qGHTGWIS9Y4HNGAEqAQBwys2bNUseOHXXUUUfpxhtv1LJly/ICsiEgh9jhh0sXXRQ08OjTx/eMACSJQBwyM2fOzPtv1apVNXr0aB122GEaNGiQNmzY8FdA3rRpU94PhIxdZ4pV2HrqKWnCBN8zApAEAnHItGvXTkceeaS23XZbXXLJJXkB99prr80LyI888kje77Es+R//+IeWLl3qe7pIJSt3ee+9wetevaRVq3jGQJbjsFYIvfPOO+rZs6ceeOAB7bjjjurfv7+eeeYZLV68WAceeGBeILYM2X5fjRo18oI1B7pCwlY6rASmVdoyffsmBmcAWYdAHFK9evXSRx99pKlTp+a9njZtmkaNGqWnnnoqLxNu1apV3pL1wbESlwiPb791jTys0pZtR3z4obTffr5nBaAQLE2HTGzvt2/fvlqwYEHeErQd0mrTpo0GDx6sBg0a5C1b268de+yx+qf1t0W4WBvLG290Yzsxb3eL7Y4xgKxEIA6Z2BLzLrvsov/973+aOHGivt58v/S5557TjBkzNGTIkLwfe+21V16ARghdcYXUqpUbf/opy9NAFmNpOsRZ8a+//qqTTz5ZFStW1NNPP513palx48Z6+OGH87Ji2yeuZLWmEU4ff+yaetjfh622siP10k47+Z4VgC2QEYdQrJpWnTp1NHDgQC1atEhdu3bVN998k7dkXbNmzbxgTRAOuf33dyenzZo1rkMT19aArEMgDrn99ttPZ5xxhiZMmKAOHTpo3333/aviFiLg1lulHXZw45dflp54wveMAGyBpekIsMNakydP1nbbbafdrJk8ouXZZ6XTT3fjevWs/Jq07ba+ZwVgMwIxEAUWiC0gm65dpeHDfc8IwGYE4jD680+pQgXfs0A2mT9f2mOPoNLWO+9Ihx7qe1YA2CMOoSVLpLfecmPujiKmUSO3Xxxjd4vXr+f5AFmAw1ph869/SUcfLR11lPTzz75ng2xiPYrtJLWxfeI77vA9IwAsTYfM7Nlu+dGWpmvUkL75xh3OAWI++0xq3dr9Halc2b3edVeeD+ARGXGYXHed+wZrrHQlQRhbsmpbV17pxr//LnXvzt1iwDMOa4XFRx+5Kkpm++2luXNdVgxsyYp7WFOI775zrx94QLroIp4T4AkZcRhYtaSrrw5eW8F/gjAKY+Uuhw4NXv/jH9LixTwvwBMCcRi89JL09ttu3Ly5dMklvmeEbHfccdK557rxsmWuSQQAL1iaznW2J7z33tIXX7jX48ZJHTr4nhVygWXBVmnNArF55RXp2GN9zwqIHDLiXPfww0EQtqspZ57pe0bIFXaY7667gteXXur2jwFkFIE4l61bJ91wQ/D6zjut9ZLPGSHX2CGtww93Yzu89e9/+54REDkE4lw2aJA0b54bn3hi8A0VSJa9cRs2zN0pNv/7n7tbDCBjCMS5yvb1/vvf4Jvp7bf7nhFylRX0uPba4MyBlb+M3UcHkHYE4lx1223S8uVu3Lmz1LKl7xkhl/3f/7mDW7E76fHXmwCkFaemc5EtR9s1JSvaX6WKK23ZpInvWSHXWUem2PbG1ltLX30l7bCD71kBoUdGnKuNHWKdc/r0IQgjNQ47LLiDbu0Se/fmyQIZQEacaz7/XNprL1dNq1Yt19ihdm3fs0KYzh7YEnWs0tYzz0inneZ7VkCokRHnmmuuCYr02wEbgjBSadttpXvvDV5fdpm0ciXPGEgjMuJcYmUs27ULGr3b3nC1ar5nhbCxN3p2He7ll91rW6IeOND3rIDQIhDn0jfHgw5yJ1rN6NHShRf6nhXCyop77LmntHatux734YeuchuAlGNpOlc89VQQhK2F3fnn+54RwmznnaWbbw7eBNrd4g0bfM8KCCUCcS6wb4CxggvGindUqOBzRoiCyy93BwPNjBlS//6+ZwSEEkvTucCKK/TsGVwxeestakojM6ZMkQ44wGXFdh5h5kyXLQNIGTLibPfbb8ESoaGxAzJpv/2C+8S2X9yjR3BqH0BKEIiz3T33SIsWubH1GbbsBMik//zHndKP9Sx+/HGeP5BCLE1nMyuq8Le/uazY9oS//FLaZRffs0IUPfdcUNijbl1p1izusAMpQkaczW65xQVhY6dWCcLwpX176Ywz3PiXX6SrruJrAaQIGXG2stKVVmrwjz+k6tWluXOl+vV9zwpR9tNP0u67uzrUsQIzdngQQJmQEWer665zQdhceSVBGP5ZJyZrvxljqzSx5iMASo2MOBtNnepOq8b24yw7trZ0gG9//im1bRsUl7nxRummm3zPCshpZMTZxq6GXH11YstDgjCyhR0aHD5cqljRvbYM2Q5uASg1AnG2seshb7zhxnZi2pb/gGzSqpX0j3+48e+/S927Sxs3+p4VkLNYms4m9s1sn32kzz5zr+2+ZseOvmcF5GfFPazm+bffutcjR0oXX8yTAkqBQJxNxo6VOnd24zZtpI8/lsqzaIEs9dpr0rHHunGtWm6Jevvtfc8KyDl8l88W69ZJ11+fWMqSIIxsdswxUqdObrx8udSvn+8ZATmJQJwthgyRfvzRjY8/XjrySN8zApIrwVq7ths/9pj08ss8NaCEWJrOBpZN2MGsX391XZU++SRoPwdku9GjpYsucuOddpK++MIVoQGQFDLibHDHHS4IG1vqIwgjl1x4odSunRt//31itzAAxSIjzoaygc2auT3iypWlr792WQWQS2bPdtearNKW3TW2ojR77+17VkBOICP2zSoTWRA2l11GEEZusoYkVpY1Vn3L7r/bfwEUi4zYJ2tr2LKluz9cs6YrZVmnjtcpAaVmxT0sC/7qK/d64ECpd28eKFAMMmKfrrkmqEj0f/9HEEZus60VK38Zc+210rx5PmcE5AQCsS/vvSdNmBB0tenTx9tUgJQ55JCgLKv10iYjBopFIPbV2CG+sbqdMt1qKy9TAVLu9tuDClvPPSc98wwPGSgCe8Q+2DemM85wY2u0brWlY91sgDB44gnp739344YN3b7xNtv4nhWQlciIM+2PP9zecHz2QBBG2Jx9tnTCCW68YIHbLwZQIDLiTLPDLNY2Lraf9s47rpoWEDZW3GPPPaU1a9zf8cmTpQMP9D0rIOsQiDNp9WqpeXNp4UL3+v33pYMPzugUgIy6++6gd7Fd1Zs2TapUiS8CEIel6Uy6994gCJ9+OkEY4de3r+uxbT7/3DWJAJCAjDhTliyRmjaVVq1yJQCtMP5uu2Xs0wPeWLnLAw5wd+arVnV/963JCYA8ZMSZ8p//uCBsLr6YIIzo2Hff4J68lXPt0cNd4QOQh4w4E777Ttp1V2nDBqlaNWnuXHelA4gKK+6xxx5Bpa2HH5bOO8/3rICsQEacCddf74KwueIKgjCip0YNafDg4HW/ftLSpT5nBGQNMuJ0mz5datPGja2hgzV2sAYPQBR16CA99ZQbd+kijRrle0aAd2TE6WbNHGJuuIEgjGizjkyxClujR0tvveV7RoB3ZMTp9Npr0rHHuvHOO7syf1WqpPVTAllv6FCpZ8+gj/GMGe40NRBRZMTpYlc1rr468dQ0QRhwleUOOsg9idmzpdtu46kg0siI0+XRR4NToVbQwO5Slud9D5DH7hLbvwurvW6Vtj791J2qBiKIyJAO69dL110XvL7jDoIwEK9Fi6AVqN0osCzZVpGACCIQp8P997uC9+aYY9wPAPmv9TVr5sbvvSc98ABPCJHE0nSqrVjhyvfF7khakfvWrVP+aYBQmDRJOvpoN65Vyx1orF/f96yAjCIjTrW77gqC8LnnEoSBohx1lHT++W68fLl0+eU8L0QOGXEqWWcly4bXrnUHUL7+2l1bAlC4X36Rdt89eAP74ovSCSfwxBAZZMSpdNNNLggbuydJEAaKV7eu61scY9W2OLiFCCEjTpVZs9xJ0D//lLbe2pWytG8wAIpn3ZhOO006+WSpa1f378jahQIRUNH3BELj2mvdNw9jhTwIwkDyypWTnn46aI9IEEaEkBGnwuTJUtu2btyggTRnjlS9eko+NIDCzZkzR1tttZV22GEHHhNyFnvEZWXv4ONLWdo+MUEYSLsPP/xQu+66q6677jr9+uuvPHHkLDLispowQWrf3o133dWV7qvIij+QCV26dNFzzz2nE088UQ8//DAPHTmJQFwWVid3r72kL790r22P6/TTU/OVAVCoTZs2qVy5clq6dKnat2+vyZMna+LEiTrppJN4asg5LE2XxZgxQRC2bjJ26hNA2lkQnjZtmnr06KEpU6bo+OOPV0VWopCjyIhLa80a10v1p5/c63fflQ45JHVfGQAFsv3gESNG6I477lD16tXVt29fnX766fqbFdMBchCbmaU1cGAQhE89lSAMpNgnn3yipk2bqmbNmnlL0Rs3btTbb7+tG264QR999JE6duyYlxG3adNG1apVS1iyBnIJS9OlYaX4br998xMsT2NzIMUee+wxHXrooXrRyl1KWrBgga688kqdcMIJWrlyZV5G/L///U+HHHLIX0HYEISRi1iaLo0rr5TuuceNL75YGjkytV8VANp5553VsGFDnXnmmRo+fLjmz5+vSy65RBdeeKFatGjx154wWTByHYG4pKzPsF1T+v13qWpVV7yjUaO0fHGAKHvjjTd09OYWiZYdX3XVVTrwwANVp04d31MDUoql6ZL6179cEDbWso0gDKTFkUceqYttxUnSBRdckHc1iSCMMCIjLokZM6R99nHVtGrXdo0drJk5gLRYvXq1mjRpkpcJDxkyRDvuuCNPGqFDRlwS//d/QVH6664jCANpZteT/vOf/+SdlqaMJcKKjDhZb7whHXWUG9u78q+/lqpUSd9XBsBffvjhB7JhhBaBOBnWpHz//aVp09zrhx6Szj8/vV8ZAEAkRLagh60wL1ggTZ8uLV7szl9VrizVqye1bi01bOhapOYZNy4IwlZb+rzzfE4dQHy9d0pbIsdVjFpi+/rr0ogRriLlokWF/97tt7crE1LXCzfo6GuvDzbTrZCHFfEA4P8f9Lffutajo0a564RADorE0vSff0rDh7saHHPnlvzPN9McXaF71O2Iuaow6dW4VBmAF/ZtywrpXHaZW866/nrpllv4YiAnhT4QW72NLl2k999P/Hm7dbTvvlKbNpLVirdlafv3bDeSbBV66lRp+fLEP9N2r1UaPW5rNW+e0f8FAAWZOdNdJ9ywQapUSfr0U2mPPXhWyDmhDsRjx0rduknr1gU/Zwefe/Z0fRqK2lqyracJE6QhQ6RJk4Kft9Uvy645qwVkgRtukP7zHzdu21Z65x22jpBzQhuILYD26hW8tqzXtpEOO6zkH8v+bVtWbdtRMYMHu4AOwCN7l92qlVv6MvffL3XvzpcEOSWUgdgy4c6dg9ddu0r9+1txgNJ/zNWrpX793EGvGG4xAVl2x79mTemrr6QGDXzPCkha6I7/2htjW46OueYaadiwsgVhY3/ePo59vBh74x17Iw7AkyOPtGLUbrxihasBD+SQUGXEdjr68MODg1mWCVvwTOUhZ3taFoBjmbFtS739tlShQuo+B4ASWrJE2n1391/z/PPSSSfxGJETQpUR2yGqWBC2PWFbjk71TSP7ePZxmzZ1r+3z2ecF4NF22wU9wo0d4PjtN58zAqIXiO1uf/y/QzuYVdRy9MyZM3XWWWepadOm2mqrrbTddtvpsMMO08SJE4v9XPZx7ePHWGC2zw/Ao06dgr3iH3+UbryRLwdyQmgCsVXMihXrsH+LxZ2OtiLyq1atyutzOmDAAN1g1yBk15pO1fAkUlxbAretKWP7xPFXnAB4YMtVdmo6VmHr3ntdDVsgy4Vmj/iss6Tx4934qaekM84o+cf4888/1aZNG61bt06zZs0q9vfb5+nQwY3tv1aSGoBnt90mXXutG1vFng8/pB41slooMmJ7K2G1o2MVs6xYR2lUqFBBjRs31vItS2oVon37oCXxe++V7nMCSLF//ENq0cKNrUzeoEE8YmS1UARi66IUa+BgZStL0oxl9erVWrJkib755hv1799fL730ko6K7TMVwz6PveE2P//s5gHAMyt3adtLsZOatu1ke8ZAlgpFII7fBooFxmRdeeWVqlu3rpo1a6Z//OMfOv3003Xfffcl/efjP1+sUyIAzw46SLr00qAaj5XZC8cuHEIoFIHY+gnH2LWlkrj88sv12muvacyYMTrhhBPy9ol/t+4PSYr/fL/8UrLPDSDNe8WxClt2r9gOdQBZKBSBOD5uWhelkthtt9109NFHq3Pnznr++ef122+/6ZRTTlGyZ9iqVAnG69eX7HMDSCMrdxm/P9y7d/6WakAWCEUgjg++JUhmC9ShQwdNmTJFs2fPTur3xwff+KAMIAvY9YlTTgkOcsTXqAWyRCgCcb16wdj6CZfF2rVr8/67wmrWJiH+89WtW7bPDSDF7MCWnfmIVfexe8aTJ/OYkVVCEYhbty75ganF8RvLm23YsEEPPfSQqlWrpj2SbDAe//lKelAMQAY0aSLdemvw2rrClHXpDEihUBT0sP8DO5NhV5jsXq8dmiruCpOdjl65cmVeWcsddthBP//8sx555JG8Qh533323rrjiimI/7x9/uCzYtp3q15cWLkzd/xOAFHeEOfBAaepU9/o//5Guu45HjKxQPiyrT4ce6sYWFCdMKP7PdOzYUeXLl9fQoUPVo0cP3XPPPWrUqJGee+65pIKwee654OzHIYeU5f8AQFpZezS7Wxxrk3bLLfQwRdYIRUZsXn1VOu44N7Z6HFZ7Ot3s81hP8rzPf8mTOubu46Vttkn/JwZQ+qpbd9/txlYs3r5RpLpFGxDVQGzdj3bdNWj8YD2Ci2v8UBb28du1c+Pmmq1Z2k3lt9la6tFD6ts3uL8IIHtYcY8997SuL+71mDFS586+Z4WIC8XStClfXopfUb7oIvdvLh3s49rHj+mne1Vem6SVK6U77pB22knq2lX6+uv0TABA6djp6aFDg9f2TWPJEp4mvApNII4dhmzbNrhW1K9f6qva2ce7/HLp22+DveHuX/aVLrkkuNBsJzJHjpR2393dY7TuLwCywwkn2CERN1661Orc+p4RIi40S9Mx1hu4VStp3Tr32u7v282FVGwD2ZOy7mq33+5eV6smzZghNW+++TfYsekBA9w7bsuO49k6+VVXSSeeyJ4U4JsV97A3yrHTlrZXnGSzFyDVQpURGwuKdjgyvtxs9+5lX6a2P28ZdywIm2HD4oKwsX1h+w3z5kl33pm4T/zOO9LJJ7t3CQ89ZJeWyzYhAKVn9w3t32iMNYjYXMwHyLTQZcQxQ4a4hisxTZtKo0ZJhx9euoNZticcW442gwdLPXsmUf/ykUeku+6SZs1K/LXGjd3aue0l16hR8kkBKPsJT/uGEGsmbstd8YU/gAwJbSA2Y8e6bDj+ja7dWLAA2r590UU/LGG1+8gW0GNXlGLL0ZYJn39+Cf/BT5zoDnJ98EHir227rZtQnz6JtToBpN9XX0l77eX+wds3hE8+kVq04Mkjo0IdiGN7xl26SO+/n/jzVoHLSlLaD2tlaA0bLIG1Q15WttJ+bNmoxQ6CjR69xXJ0Sdm7b1sSs8Acr2pV6cIL3T3HkvZyBFB6//qXK/BhDj5Yevdddw0DyJDQB+JYdTvbN77nnuCecUk00xxdcX11dbup4V+Fecrsyy/dkrUtXcfvF9s3gDPPlK6+muLVQCbYyU7LimMd1+ywpe0ZAxkSiUAcv0I8aZILypaY2sHJos5yHFJ/jrp92lNHaZLKX3KxNGJE6ic1f750771uvfu33xJ/zdbRLSAfcwwnrYF0eust6Ygj3Niq49mSdcOGPHNkRKQC8ZYWLHBL0NYkwpalbXnamjjYcnXev8FVq6RGjdxVJPtFOw2drl6Hy5a5Fm12/cm6V8Tbe2939emss4rvZgGgdOxEpu09mQ4dpHHjeJLIiEgH4hLXprV9pOuvT/8ymV1vsmXrLdfRrWKXVQKybxix/qoAUsOKe+y2W1Bpy85x2JVDIM0IxMX5/nt3eMrWtW292mrUxipopXtj+9ln3UnrKVMSf61OHemyy9yP7bZL/1yAqLAzG506BVcM7SwH1wuRZhwNLI5loaef7sa2qfzEE8oIOxVmh7Y++kh6803p+OMT37nffLNreN67t3uzAKDszj1XOvZYN7atqBtu4Kki7ciIk2F3n2INh/fZx20s+2id9tlnbsn6scdcxhwftM8+2+0j234ygNKzO4x2l9i2iewWw8cfc4MBaUVGnAy7W7jvvm5sF/5jlXgyzcpjWpUS+0ZhrRa32sr9vAVlC872JsGaMlsFErb+gdKxragbb3Rj25Ky6nd//MHTRNoQiJNh2a+1XIqx60Y+7bijm8OPP0r//nfiPvGrr7ri9fvv7059xmfOAJJjHZlatgzefA8cyJND2rA0nSxrbWj7xdZhyZar7ETzzjsrK6xZIz34oPS//0nffZf/3b2d/L7gAlefE0By7HzGQQe51SVbfZo5030PAFKMjDhZdlLaTinHlqsGDVLWsG8SVq/aKgM9/rhboo6xZewePdw3ECto/+uvPmcK5I4DDgg6u9ibXesiw5YP0oCMuCTsfqFdabBDHFtv7apiWRWebGPfLKyEmF19sj6r8ez+sfVztM5P9v8CoHBWzMf6Flv1H2O3JuxgJJBCZMQlYXuxsbZLVnUrVoUnG/e0jz5aeu01d8K7Y8egiL01Vu7f3/WFtOXqL77wPVsge9kb7fvuC15blzSrggekEBlxSdkF/z33dGMLZrYcnLJOEGlkzZStQpg1ZbaMPt5JJ7mrT4ceSk1roCCnnSY995wb24qS1YYHUoRAXBp2RchOJxurfmXNjXOFFda2d/j2Y8v94gMPdAHZ/n9oAwcErLjHHnsEjVmsVWKstgBQRixNl0Y2XWUqKWtaYVW5rFSnNZiw6lwxH34onXGG2xMbOdJ1wgDgzlPYYceY7t3dTQogBciIS8NOTdu746+/Du4Z5mpFK+uF/OST7mDX558n/prV1rY3HdabtWZNXzMEsoPdybfrTLHa75loAoNIICMu1VMrn9tZcbxKlaTzzpNmzJBeeklq1y74Naut/X//57IBW7L+6SefMwX8srMg1pM8dibkP/9xZ0SAMiIjLi07fWwByk5Q2h1jW+q1DDIMrLbunXdKTz+deG/SgradGrcCIbZ8DUSRvSm1mu/miCPcVUEftecRGmTEpRW7j2tsr2joUIWGlcccP94tvdv/Y5UqwTK2nbq2ZXk7RTp5su+ZAplndahjFbasM9qYMXwVUCZkxGVhBT3sH6TtHdkhKKv9XLWqQseWqK3W7pAh0ooVib9mJ0ctQ7ArUJy0RlS8/LJ0wgluXLu2NGuW+x4AlAIZcVk0aiSddVZwLcg6IIWRLbn/97/uCofVs95hh+DXrBPVqae6AvlW75qTpIgC6w9+zjlubNcAr7jC94yQw8iIU1EY3u7fxtoUfvpp+PeLLNg++qjbJ7MCJ/EsSFv5TFvStjKgQFgtWiTttpu0fLl7bbUFjjnG96yQgwjEqWBXGuwOrrGDG0ceqchc43rhBXewa8sezXbdyQrmW0nAsBxiA7Zk9+2tX3Gs05ldAaTLGUqIpelUsAwwDFeZSsr2hE85xVUZev/9xApjtpd8221uD92KH8yZ43OmQHpcdJF02GFBpzO7WwyUEBlxKvzxh6s7bXuotixtp42bN1ckffWV20ceO9adso6x52JVu+xgl53KBsLCDmrttZfbsqlYUZo+3Z2ZAJJERpwK9o8v1qvY7t3aCeOosvvFDzwgffed9M9/BvvE9lyeesr1eLW7l1Y8hN6uCAPbJ77mmuBNuZ2PsG0bIElkxKlihT3sFLU1ELc7xna1qVatlH34nGUHWaxTjS3Z2zWoeHa4zTJk6+9qxUKAXGV12S0rjpW9HTzYnZEAkkBGnCrbbitdeGFQdcuyQrg3I1df7TJkKw+4yy7BU/nsM6lTJ6lZM7eKYM8NyEVW9Ca+NaJlyJSERZLIiFPJ3g3bMpXZcUdp7ly3bI2ALdlZX1drMmFXv+JZYQRb4rcfFEdALrrkkuBNuJ2JsO0YoBgE4lSzClMvvujG48ZJHTqk/FOEgu0P22lru/pkV6Di2fUPO41qRRLsEByQK6y4h52TWLzYvbY3nVbwBigCgTjVXn89uNTftm3++7XIz+5e2klrKxJih13ir0fZ/rEd+mrdmieH3GB/j62jmbFzI1b0huI2KAKBOB2Znl1dmDkz6GS0334p/zShZLW67VDX8OH594uPPtrtNR91VPgrlyH3vwdYHepXXnGv+/aNVn0BlBiBON3Vds49V3rkkbR8mlAv71k3qwEDXA3veJYZ20nrM89k/x3Zyw4n7rmntHate+No5yF4Q45CEIjTwf7xNWkiLVnigsX33yc2SkDyz9FazNmytVUtimd7x1de6U6qb7UVTxTZx84/2CqO2XtvacoU3jyiQFxfSgc7bHTppW5se57WPhClf452Gv3JJ6U2bYJf+/ZbqVcvdzrdygouXcoTRvaVvrW78saawbA8jUKQEafLggWuzrKVeaxTx+1/krmVfe/NGrFbphHbf4uxZ2vbAfbNz4IzkA3sjIh1Z7O/u/Z39IsvpJ139j0rZBky4nRp2FDq2NGNLVt7+OG0farIsL0262xlTdk/+cTtv1eo4H7NKprZnrJ1wDn/fFcsBPDN6qrHyt/a31GrtkVpV2yBjDidpk2T9t3Xje1uoZ2k5sRv6g/F3HOPK6Jge8rx7OSqHew6/HCeO/xZuVLaY4+g0tZjj0l//ztfEfyFQJxu1iLNClcYW0499ti0f8pIsoNx993nfmy5X2ynVe3QzGmnBRk0kEnPPiudfrob16vnOjZZWVyApekMuPzyYNy/P3/p0mW77aSbbpJ++MHVrY7fJ7bTqlbhzFYl7I7yunV8HZBZ9iYwFoit6lbsNDVARpwBf/7pmhrYFSZjVXYsICC97LS6nbS2g10zZiT+2vbbuyILPXrQIQuZY0vT9m9/1Sr3+p13pEMP5SsADmulnS2F9ukTvI5yr+JMsvvbdpjLDnXZ4S475BWzaJF07bVS48bSP/7hWlYC6Wa1BP773+C19S229omIPPaIM2HFCldz9rff3N1Y+8ZvnYaQWVOnugzZOuLEN263XshWG9hqWtuhGiCdK2RWgz7Weezmm6V//YvnHXFcX8qEmjVdNyFjJ3ttnxKZZyfYbbnaCoRYoRDrIWvsrveDD7qShNYph0YdSOcKmf37j7VHvfVW9/cRkUZGnClWorF5c3eH0Jao7NqNZWLwx5aoBw2SBg+Wli9P/LWDD3ZXn045xXWBAlLp//7P9eQ2dr3OCtVwtTGyCMSZPjlp/UkNdwmzhx2esUYddh95y/3i3XZzS9a2dB3LoIGysuIeLVq4N+TG7sHHVs0QOQTiTHrrLemII4KKOx9+yLvgbGJL1PYGyfaRY20s4yul2VW07t2lbbbxNUOEyauvSscd58Z2p9juFtsdY0QOgTiTbFna2vhZAXgzebJ00EEZnQKS/Dq9+KILyHbFJJ4FYbv2ZNefGjTgcaJsOnUK2qTaqgulcCOJQJxp1tbPWveZs8+Wnngi41NACdiqhQVkq4wUXyO4cmWpc2d3/WnXXXmkKB0r7mF3i60Ht7GrdrEsGZFBIM40uzdovYrtH6CdoLR2fvYa2c1Otlpf5Icekn7/Pfh5O2Bje/9WKemAA3zOELlq1Cjp4ovd2DozWYcmOrVFCsdBM80O/FgHltidQquNjOxnWe+IEa5CmgXd2D6xZcnPPONa3dnpV1vSprsOSqJLF/d3x9jhrX//m+cXMWTEPlg2bFWdLLOqVUuaN0+qUcPLVFCGjjrDhrn64QsXJv6anYa1k9bnnMMVNSS/4tKqlfueYCtl06e714gEMmIf7GSkHcwwdn/VljuRWywjtmBrGYxdPbFrTjG2tHjBBa43sgVqq6gGFLfict11wUqZlb+0/yISyIh9sUYEe+/txlbow64uUDgid1nJzIkTXZGGDz5I/DW7mmLbEVZznOspKOr8iH1PsO8FxorNXHYZzysCCMQ+HXWU9MYbbvz889JJJ3mdDlLESmTaSWsLzPGqVnUn5u2ktWXLwJasd7n1MDdbby199ZWrxIdQY2k6W3oV33uvz5kglQ45RJowwS1RW+CN1RW2Psj33y/tsou7ujZtGs8diawtYteuQcW33r15QhFARux7OdO+KVsdavP55+6gD8LFymbaGy073LXlfrG1Z7RT2MccU+oqa5s2bdLq1atVgwN/4bBsmbtbbLXQjZ3KtytyCC0yYq9Pv7yr0BRDVhxO1gLT7iD/+KPrR7v99sGv2daEFXCwimtWXvOPP0r84SdNmqSLL75YHTt21G233aaVdqIbucvOFMR/L7B9Yr6moUYgzoY7hNYm0Vh5u19+8T0jpPMb7DXXuLvIlh03axb8mpU9Pfdcd3DP7pZbU4Aks+GvvvpKS5cu1W+//aaRI0dql1120YABA/THH3/k/TpyUMeO0vHHu/FPP0nXX+97RkgjlqazgR3euftuN77lFv7RRYVdT7HSmXbSesqUxF+rU0e69lqpX7+kl6wtGM+bN09Dhw7VhAkT9Pzzz6tNmzbpmTvSz67G2VaVvSmzvwNWbtWaxSB0yIizgS09xa4uWW9cu8aA8LPCDWeeKX30ketHG8uAzNKlrkJXEndJN9pZg7zYXUd77723brrpJq1bt053xPrdIjdZucubb3Zj+7tgd4utQxhCh0CcDXbaSTrjDDf++WfpySd9zwiZZNlOu3bSSy+5++VW7MUyYuvyFDtxXYTym9/ExQJygwYNtMMOO+jXX3/NW65Gjt+siNUbsL8bViAGoUMgzsarTPaPjb29aLKyhnZWYM4cd++4BGIB+ZVXXtHChQvVsGFDTlLnOnsjNnx4sGJ2001uyRqhQiDOFgcfLO27rxt/8okrCoFoH+xKotLaloex3njjDfXs2VO1atXKO0mNENhvv+A+8dq1bqWEN+qhQiDOpuXJLbNioBC2DG0/ym0+yGVZsF1fOvroo1W5cmUNHjxYh8c6+iD32SFOuwZnXnlFevxx3zNCCnFqOptY5xU7oLFggQvMc+dKTZv6nhWylJ2SfuKJJzRkyBB9+eWX2mmnnXT55ZfrzDPPzNsjRshYtbb27d24bl1Xk7p2bd+zQgqQEWeTypWlXr3c2Jae6FWMAnz22Wfq1atX3n3hq666Sk2aNNFjjz2m2bNnq0+fPnl7wwihU08NDnVavYGrrvI9I6QIGXG2WbLE9Sq2usRW9N3KI8aa0CPyli1bpmbNmuX998QTT9S//vUv7bfffn8tUf/555+qYNeiEE5W3GOPPYJKW2+/HTSJQM4iI842220nde4cFH0fPdr3jJBFKlasqE6dOuWNP//8c7355pv65JNP/rqmVGQQtjvJHPLJbbblcNttwWu7W0zdgZxHRpyNvvxS2nNPN7Y94tmzXfEHYDO7I3zzzTdr+PDheVeULrzwQp100kk68MADVbWga09Ww9rqFz/1lFvStL1G+l/nJrsv3ratq7RlbrzRXWtCziIQZytrBPDqq25sZRBjhzSAOFZP+s4779Tdd9+tFStWaPHixapd0AEeq8hkBwFtadNY169//lM6/3ypShWeaa6xTm3WKMTeYNnZEiv2sdtuvmeFUmJpOltxlQlJLlVfe+21eSeo7QpTgUHYlqOt81P8r9kqi/W9tapuVgpzxQqedy5p2dK9kYrdtuje3WXKyElkxNnK/lHZ8rRdUTDTp0v77ON7VshlFpDtDqoF3rfeSvw1Oxh46aXuDSCnrnODFfewgBzrZz5ypEQRl5xERpwrvYoHDPA5G4SBnay2xhLWYMIaTVjDiVhnJzsYeNddLkO2b+ZffeV7tihOtWrS/fcndnFbtIjnloPIiLPZ6tXuKtOyZW4f6IcfpPr1fc8KYWI1rf/3P2nMmPynb+1cgh3ssvKryF62z2/1yc0550iPPup7RighMuJsVr262/uJ7QMNHep7Rgib5s2lYcOk77+XrrlGqlkz+LXnnnOncw89VJo4kT3IbHXPPcH+/2OPSS+/7HtGKCEy4mxnBT1sudDugFpZOzt0U8KuPEDSrFDEiBGu1nnshHWMFZKwA0LnnutWaJA9HnxQ6tLFje37xRdfuDfyyAlkxNnOCr2fdVZQ1s7e8QLpYlXcrrxS+vZbV0xm990T77fbN3u723733W5fGdnhggukI45wY1vduPlm3zNCCZAR5wI7WHPggW5spyTtzmDskA2Q7tP7L7wg3Xln/tactozds6c7VLj99nwdfLMradbP2vb6rQDQ1KnS3nv7nhWSQEacCw44QDrooOAiv516BTLBTu+fcor07rvS++8nFpaxu8dWbnHHHd1ZBjv4BX+sSMv117uxbWVZ+Uv7L7IegTgXC3xYqUIg0+z0tFV5syXqiy6SKlVyP28Z2PDh0q67Sh06SFOm8LXxxU65216+sa/DkCF8LXIAS9O5wkrZ2d7cvHluWfrrr92JV8AXO8xl99vtLuuW+8Xt2rmgYPeW2UbJLNtCsJPupkYN98bJrkEia5ER54qKFaXevYMKSQMH+p4Ros46AdnesZ3kv/32xDvuVrnrxBPdHuUjj7ha18iMQw5xy9LGunLFvm8ga5ER5xIr7GGnqNescVcT7GpTrVq+ZwU41kPbCktYhS47OBSvSRN3GtuqdnGtJv2WL3cn3n/+2b1++mnp9NP5m5qlyIhzybbbShdeGFTdeuAB3zMCAna//ZJL3FKofeO3Q4YxljXb6WoLyNa2z67iIX3sDXp8WdzLLnN3xJGVyIhzje0Nx9qd2Tc1K/huy9ZAtrEtFDttbcvXdgVqyzrJduDriivc2Qek5/mffLL04ovuda9e0n338aSzEIE4F9k/rtg3tnHj3ElVIJvZtTuraW11kO3gYfz1qLPPdge76C6Welbcw7q42XaWHZqbPDmoSYCswdJ0LuIqE3KNFaKxxhK2gtOvX7BPbAVDHn/cNbk/9ljp9dddJofUsHKXt9zixvZc7RAXB+eyDhlxLrJ/UFZBx+rJmo8/lvbbz/esgOT9+qtrYmL7mFvuF1tQtgzZ2jSy7VJ2tgJh+/XW09zYCferr+ZvaxYhEOcqO6hlB2OMFeG3KyJALja3t0zZlq1jDe5jbO/YTlrbAcWttvI1w3CYNk3af3+3AmGH6uxN/N/+5ntW2IxAnMvfwOyw1pIlLmuwvSC71wnkIivFaCet77jDBY14220n9enj6lrXqeNrhrnPtgRiVfmOOUZ65RWKrWQJ9ohzlZ06vfTSYOmJUnbIZdakwLqMWVnGSZPcfnGMvdn817/cG087H/HDDz5nmrtsr9ieoXntNXdwDlmBjDiXLVzoCu7b4QtrDG7lL1nCQ1h8+qm7+vTkk4nNCyxon3OO641sZyWQvOefd008jPU3/+orVhmyABlxLmvQQPr734PDL1bVCAgLK49pWZt1dbKCFLYKZCwo29/1vfZyZTStnCYnrZO/+hi77miH5OzNDLwjI851tp+2775ubCXtZs5k3wfhZEvUVpDCfixdmvhrdhDJTlqfdprLmFG4BQvc94pYpS1rq2pNOuANgTgMDjvMVTAyL78sHXec7xkB6WPlXUeNku6+O/9+sXUk+8c/pM6d3elgFMw6ZvXoEfQxnjGD5+URS9NhOQ0ZQ69ihJ0VA7GOQnPnumt7tkQdY8vY3bu7Qha33eaaHyA/K+xh/aWNNeiwZwVvyIjDwPbMmjVzV5iMFd23pScgCmx/+NVX3cGuN95I/DXrx2uB2U5bW+cyBOwusZUVtVsXlSq5w3F77MET8oCMOAxsT8zuWcbQqxhRYjWUbTvGrj1ZlTm7BmU1rGP9eG0J24qDdOni3qTCadHC7asbu3lhb1is4Acyjow4LFascO/47RuPnS61q0wUP0BU2bK1BeDRo6X16xN/za7vWAA65BBfs8uuwkB2Bcyelxk+XOra1fesIoeMOCxq1nRN12P/uEaM8D0jwB/bqrFa1naY67rrXH/emIkTpUMPldq2lZ57LtpZoL1pt4NbMfYG5eeffc4oksiIw8Rq9dqpUdszs3KX333n9n6AqFu1Sho5UrrnHmn+/MRfs/7edp/2vPOkKlUUSRdcID30kBt37Og6YiFjCMRhY/co7V2+sWIIVoEIQLAX+thj7mCX3bmP17ChO9Rle6XbbBO9O9r2hiR2P/vFF6UTTvA9q8ggEIfN228Hl/OtNeJHH1HgA9iSrRpZsLGA/M47ib9mQdju2Pbt66rXRYVlxJYZGyuda29UYn2jkVYE4jB+g7F+rnYVwUyeLB10kO9ZAdnrww9dQH722cRSmZUru8IgViBk110Vevb/bl2Z7PS5sf/vu+7yPatI4LBWGK9y2PJaDAU+gKIdeKBrwWgNEKzHtwVg8/vvbl/Z7uSfcYZbXQr79w474BbbJ+/fX/rkE9+zigQy4jCy6xq2tLRokbtj/O23QfszAMXXYra7+BaUYvWY48vJXn212z+1wBVGt94qXX+9G1sde1sxoH53WpERh5G9o7Um6rGqW1YkH0By7NDW7bdLP/7olqzj94ltP/mkk9zd27Fj3eGvsLET5Hvu6cZTp/L9IwPIiMNq8WKpcWO3vGZ3KK3Ah5X7A1DyFSaraW37pbNmJf6a/Ru74gq3pB2mf192tsTuWRs7sGUVyVhVSxsy4rCqV8/dizRW+H7MGN8zAnJ3hemii9wpYjvQFX/40d7gWtMVC1I33ODeAIeBNYS49NKg25X1g6bnc9qQEYfZZ58FnWms0Ie9m4/V4AVQeu+955atrUpXPGu9eOGF7sTx3/6W20/Y3sDbQbVYpa3x46Uzz/Q9q1Diu3KY2T7WkUcG7eFeesn3jIBwsDrVEya4DkYWeCtWdD+/bp0rGWk9fs8+W5o2TTnLtrQGDQpeW+tJq2mPlCMQhx1XmYD0sUNN1ljCysleeWWwT2z1q8eNc6eOjzrKtWnMxaVdy4BPPtmNFy6Urr3W94xCiaXpsLNvCFaMINZdxZarW7b0PSsgnJYtc9eeBgzIv1+8996uqYK1aYxl0LnATo9bn2LbK7YrW++/T5GgFCMjDjvbE7ZSfTH2DQJAemy7rcsarevTsGGuC1SMVbs791x3XsOuFK5ZkxtfBTuIdsstbmxZfbdu4by25REZcRRYj2LrVWz7O3YC1E561q3re1ZA+Nk9fjtpfccd0pQpib9m/cJt37VXL2m77ZTV/vjDVSCL7Xn/97/SNdf4nlVokBFHge1b2T3H2J3I+P6jANLHKlLZPquVx3zzTen444Nfs05HN93kquD16SN9/332fiVsKX348ODWxb//HWx3oczIiKPClsqaNnV7xvXru3/0Ue29Cvhk5zTs6pP1/LWMOT5o20lr20e2/eRsZAfSrKezsUNor70W3lKfGURGHBX2rtsK1xu7F/jkk75nBET3WuHDD0vffOMy4a22cj9vQdl6Je+zj3TccdIbb2TfSeubbw4qbFmXJvv/QJmREUeJnXa0+4/G/rHbfg/vZgG/bIl68GB3Z3fJksRfs+tPliHbm+hsabxgfZyt3raxvW3rWpXte9xZjkAcJfbuev/9XSF38/bbrpsMAP/sFLXdSb77bncvOZ5V6bJqXRdcIFWrJu86dgxW1aygic0bpUYgjppHHw1qUJ9+uuvDCiC7Tig/9ZQ7ab1lP2CrIW/L2dZdza5K+WLbW7vtFlTasmXqWBU/lBiBOGqsG9POO7ueq7YsbScf7RAXgOxbwbIAZwH59dcTf806Itl9Xms4YR2gfLB70rHGEHY32g6hWa1tlBiHtaKmcmV3bzH2Dz2+liyA7GFvlI8+2p1MtvMcthwcuz5kVa7693dvom252mpeZ1rXrkGrRKtlf+utmZ9DSJARR/VwiBX4sAL1W28tzZ8vbbON71kBKM6337o95FGj3L/feHaAyg52HXpo5g5hWp9iu2pllbYqVXJL6VZ/GyVCRhxFVtGnc2c3XrWKgxZArrAM2E5YW/1n638cv0/8wgvS4Ye7XsLPPONqBqSb1aC++mo3tmBsy+WZ+LwhQ0YcVfZONvbO1f5xz56dPdcjACRfvvaBB1yRDQvO8awV4z//KZ1/fnqL91hmbnejbXnaWOW+7t3T9/lCiEAcZVY0wNqzGXsHfdppvmcEoDQsG33iCVex6/PPE3/NKulZO1Q7WFWzZnqer5XvjJ2ats9hd4sbNEjP5wohlqajzE5cxtx7r8+ZACgL25/t1EmaMcMV3GjXLvGq0f/9nztdbXvIdmMi1Y44wt0nNnalKb4POopFRhxltpdjy9OzZrnX06e7ilsAct/HH7sM2WoFxJfKtKBty9W2bG13gVN5CNQ+Xqw62PPPBxW4UCQCcdTZfk6PHm5sB7jGjPE9IwCpZHu3//uf+7dt3dfitW/vDlsddFBqPtfYscFBUKtJPXOm6/62mb0fsITc3vMvXuzKGlSu7OqUtG4tNWwYzaq7BOKos7J6dpVp2TL3TtkOfNieEoBwsSXqgQOlIUOCilgxVoPeAvKJJwZ3lUvDIu2xxwYFSK64Qhvvujvv5YgR0rvvSosWFf7Ht9/e3b6yK8p2hbosU8klEfnfRKGs80vshKMd+Bg6lIcFhJG9wf7vf92bbcuQd9gh+LX33pNOOUVq2VJ68EGXqpaGpbO2yla1qv5UeQ3tv0677rQu71zo+PFFB2Fjv26/z37/rru6b0fxnSLDiowYrqDHTju5v/F167p/qJSqA8LNgq3Vnrd9ZDvlHM+CtB3mtHvBVvSnhOZceb+63NNC72tzt7fNatVyDaXatHF9LGxZ+vffXUdIKx5m/WiWL0/8WFa8y3pKWBXNsCIQwznnHNeo3Ni9xIsu4skAUTm0acVArKa1tUqNZ1eRrMFE375u3TjJbeJu3TZp3bpgs/eoo9yHOfVUqWLFovtdTJjgVs+tzHaM5QXDh7szZmFEIEZwwvKAA9zYlqfsGkQUT00AUWaB+K67pOeeS/x5KwhiNa2tFWMRqakF0Fgpe2NZr1XjLE231Xfekbp0cVU9Y6yomAX0sCEQI2Cl8T74wI1pawZEly1V2z6ypbd2diTG3pyfcYY72LXffoUemDZ24Mr6UlijqNJavdqtkNtBr5iHHgpfZkwgRsAafVuHF3PyydLEiTwdIMp++kkaMMAdwLK69PGsaIgF5OOO05y55fKqXMb6UFxzjWvGlIpFtU2bpOuuk267zb2uVs0t2IVpz5hAjMQNGqs7PW+ee231p8P0tx1A6dgJKus/bBX47BpUnD9b7q3D17+i92fX+ysTtt+ayp2tTZvc5Y5YZmwHuN5+Ozzl8bm+hICdoujdO3htdw4BwI47W/b73XcuGlpDic2Gf37gX0HY9oRtObqoILx+/XpdffXVatiwoapVq6YDDjhAr1nP5SLYx4u1X45tZdvhrbAgI0YiK+xhBT6s0Idt7tjVJvtHCAAxdtVxwgRtvP1O7frxQ5ort3JmWWpxB7POOeccjR8/XpdffrmaN2+uBx98UFOmTNGbb76pQ6ywSBHs48fKaNtinVXnDUPRDwIx8rNjj3b80dgJSjspCQBbePWVTTru+HJ/XVGKFdQqzMcff5yXAd911136x+bvK+vWrVOLFi1Ur149TZ48udhnbJ/njTc2f/5XpWOOyf0vSwjeSyDl7M5gzKBBbu8YALYwYmSwBp3MtSLLhCtUqKBuVihks6pVq+riiy/WBx98oHmx8ylFiP88YVmeJhAjP9v/iXVNsSpbzz7LUwKQ7wCV1Y42tntlxTqK88knn2iXXXbRNttsk/Dz+++/f95/P/3002I/hvWpiO2WWWXOMCAQo2Dx/UTtlAQAxLEuSrHa0Va2sqiKWTELFy5UgwYN8v18g80/tyCJXsn2eaxEprED3Olor5xpBGIUvhHTooUb276NVd4CgM2slWFMLDAWZ+3atapiVbq2UHVzbXv79WTEfz6rUZ3rCMQo/L5AfFZsl/oBYDPrJxxj15aSYdeV7PrSltZtrgRiv56M+M/3yy+5/yUhEKNw554rbbddUHXLquwAwObmTTHWRSkZtgRty9NbWrj55+xucTLik+oC4nrOIRCjcPbu9NJL3dhOTlvFdQDYIvgm275477331uzZs7Vy5cqEn//oo4/++vVkxAffAla6cw6BGMXfFahUyY2tbp0V+gAQefVcMa081k84GR06dNCff/6p4XH3jtavX6/Ro0fn3S9u3LhxUh8n/vNZC/Vcl8Q5N0SanWb8+99da5Vff5Ueftg1CwcQaa1bl/zAlAXbs846S9dcc40WL16sZs2aacyYMfr+++/1gPVBT1L850v2oFg2o7IWkjseGfvbvvvu0syZ9CoGIs7uEdv7dLvCZPd67dBUMleY7GDWDTfcoIcffljLli1Tq1atdMstt+i4445L6vPaLpllwdaHon59219WziMQIzmHH+46dZuXX85rfQYg2s46y6plufFTT7lWxen21FO2xO3G9t9x45Tz2CNGcuKvMlkrNACRZy0PY2Ll6dNtSNznCcsuGRkxku+20qyZ9P337vWXX7plagCRtXGjtOuu0ty5Srr7Ulm8HdLuSyH4X0BGWAfuPn2C1xT4ACLPguAVVwSP4aKLpNWr0/NYVq92Hz+mX79wBGFDRozk2d0/61W8apW7Y2ydUurU4QkCEV8ssyMk778fLFfbTUcrzpfKg2HdukkjR7rX1rbYsuOwBOKQ/G8gI6xjSuwtqdWEHTGCBw9EnC2WjR5t9aLda/u2cN11LnimwqZN0rXXBkHYcoBRo8IThA0ZMUrGbtLb5oz969hhB+m774KCHwAiy0oNdO4cvLbM2Bq3Va9etuXoyy8PgrB56CHp/PMVKiF6T4GMsGrr1hDUWO3p2N0FAJFmwTG+Cq5lxq1auSXk0nj7bffn44OwffywBWFDIEbZexWnag0KQM5XxLWMNdZE6dtv3Sln66pq93+tGEdRNmxwv89+v/05+/PGPp59XPv4YcTSNErOAq/Vt/v006Bf8UEH8SQB5JkzR+rSJTjAFWMVuKxIn/2wxTVr2LB+vdvxsrKV9sMqZsVr29btQduOWFgRiFE6Y8ZIF14YlNexNokAEHea2no73HNPcM+4JJo1c1ej7LS0HQgLMwIxSsfexu64oys0a8cXbQ3JXgPAFkU/Jk1yQfm996Sffy788dSv764mWfC15ekwnYwuCoEYpffvf0s33ujG//yndOedPE0ARVqwQJr25Df6pd+tWq8qqtLuYNW94vy85eqGDaP58AjEKL3Fi6UmTVx2XLOmNH++VKMGTxRA0aZMkfbf340vu0waNCjSTywiiT/S1hn8vPPceMUKt28MAMlsIMdUTKJ3YsgRiFE2ffsm1p+2DSEAKEr8PaYKIT+JlQQCMcrGbtwfeWRwZ+Gll3iiAIpGRpyAQIzUF/gAgKKQEScgEKPsTjrJXfozdk/h8895qgAKR0acgECMsrPLflvuFQNAYciIExCIkRpWZcuuMJmHH3ZXmwCguIy4Aoe1CMRIDbs/bH3PjN0rts7gAFBcRlyR60sEYqSOXcyP1aQbMsQFZADYEhlxAgIxUsdqTZ9xhhtbQVkaQQAoCBlxAgIxUotexQCKQ0acgECM1Dr4YGm//dz4k0+kd9/lCQNIxPWlBARipFa5colZ8b338oQBJOL6UgICMVKvQ4egn9mzz7pexQAQQ0acgECM1Ktc2Z2gNps2Rb7FGYAtkBEnIBAjPbp1k6pWdeMHHpBWruRJA3DIiBMQiJEedepInTu78apV0ujRPGkADhlxAgIx0mfL+tPx74IBRBfXlxIQiJE+e+whHXecG3/3nTRxIk8bAAU9tkAgRnpxlQnAlsiIExCIkV7HHivttpsbv/22K/IBINo4rJWAQIz0siYQZMUA4nFYKwGBGOl3/vnSttu68WOPuYYQAKKLjDgBgRjpt9VWUvfubrxhgzR0KE8diDIy4gQEYmRGr15BA3ALxOvW8eSBqCIjTkAgRmY0aiSddZYb//KL9OijPHkgqsiIExCIkTlbHtqyOtQAooeMOAGBGJmz//7SQQe58eefS2+8wdMHooiMOAGBGJnFVSYAFPRIQCBGZp1xhtS4sRs//7w0Zw5fASDKGXHFzYc4I4xAjMyyf3S9ewevBw7kKwBEDRlxAgIxMu+SS9zdYmPtEZcv56sARAmHtRIQiJF5VmWrSxc3Xr1aGjmSrwIQJRzWSkAghh99+gTjQYMS/2ECCDcy4gQEYvixyy7SSSe58Y8/Ss88w1cCiAoy4gQEYvjDVSYgmsiIExCI4c9RR0ktWrjx5MnSxx/z1QCigIw4AYEY/pQrl5gVDxjAVwOIAq4vJSAQw69zz5W2286Nn3xS+uknviJA2FHQIwGBGH5Vqyb16BH84xw8mK8IEHZkxAkIxPDPAnGlSm48bJi0Zo3vGQFIJzLiBARi+NeggfT3v7vxr79KY8f6nhGAdCIjTkAgRnZeZdq40edsAGQiENuBzfKEIZ4AskPr1tJhh7nxrFnSa6/5nhGAdC9NV6jAMyYjRlahwAcQrYyYFoh5yIiRPU49Vdp5Zzd++WXpq698zwhAOpARJyAQI3vYMlV8MwgKfADhREacoNymTZs2Jf4U4NHKlVKjRtKqVe6O8bx5Up06fEmAMGneXJo7V6pdW1q6VFFHRozsss020kUXufHatdLw4b5nBCBdGTGHtfIQiJF9evd21xrMffdJGzb4nhGAdOwRc1grD4EY2edvf5Pat3fjBQuk8eN9zwhAKpERJyAQI/uvMvXvL3GUAQgPDmslIBAjO1lxj733duMpU6QPPvA9IwCpwvWlBARiZCfbI+7XL7HsJYBwICNOQCBG9urYUdp+ezd+6inphx98zwhAKpARJyAQI3tVqSL17OnG1gTCTlADyH1kxAko6IHstnix1KSJtH69VLOmNH++VKOG71kBKIvKld21xH32kaZPj/yzJCNGdqtXTzrvPDdesUIaM8b3jACUFRlxAgIxsl/fvon1p+lVDOQuu4oY+zdMZa08BGJkv1atpCOPdOM5c6QXX/Q9IwBlzYYNgTgPgRi5gatMQLhOTBtKXOYhECM3nHii1KyZG0+aJH3+ue8ZASgNMuJ8CMTIDeXLJ+4VU+ADyP1ATEach0CM3HHhhe4Kk3nkEXe1CUDuLk2zR5yHQIzcYfeHu3Z1Y7tXPGyY7xkBKCky4nwIxMgtl13mlqnNkCEuIAPIHWTE+RCIkVt23FE64ww3/vln6YknfM8IQEmQEedDIEbuX2WiVzGQO8iI8yEQI/ccdJC0335u/Mkn0rvv+p4RgGRxfSkfAjFys1fx5ZcHr/v39zkbACVBQY98CMTITR06SA0buvFzz0nffut7RgCSQUacD4EYudtGzU5QG9sjHjTI94wAJIOMOB8CMXJXt25S1apu/MAD0sqVvmcEoDhkxPkQiJG76tSROnd241WrpFGjfM8IQHG4vpQPgRi5Lb7+9MCBif/IAWQfri/lQyBGbttjD+m449z4u++kiRN9zwhAUciI8yEQI/dxlQnIHWTE+RCIkfuOPVbabTc3fucdafp03zMCUBgy4nwIxMh91gQiPiseMMDnbAAUhYw4HwIxwuH886Vtt3Xjxx5zDSEAZB+uL+VDIEY4bLWV1L27G2/Y4FokAsg+FPTIh0CM8OjVS6pY0Y2HDpXWrfM9IwBbIiPOh0CM8GjUSDrrLDdeskR69FHfMwKwJQ5r5UMgRrjEH9qiVzGQfTislQ+BGOGy//6uX7H5/HPpjTd8zwhAPDLifAjECJ9+/RKzYgDZg4w4HwIxwuf006XGjd34+eelOXN8zwhADBlxPgRihI+dnO7dO3hNgQ8ge5AR50MgRjhdcolUvbobjx4tLVvme0YADBlxPgRihJNV2brwQjdes0Z64AHfMwJgyIjzIRAjvPr0CcaDBiV+AwDgBwU98iEQI7x22UU66SQ3/vFH6ZlnfM8IACUu8yEQI9y4ygRkFzLifAjECLcjj5RatHDjyZOljz/2PSMg2jislQ+BGOFWrlz+spcA/OGwVj4EYoTfeedJ223nxuPGSfPn+54REF1kxPkQiBF+VatKPXoE78bpVQz4Q0acD4EY0WCBuFIlNx42zN0tBpB5ZMT5EIgRDQ0aSH//uxv/+qs0dqzvGQHRREacD4EY0bHloa2NG33OBogmri/lQyBGdLRuLR12mBvPmiW99prvGQHRQ0GPfAjEiG5W3L+/z5kA0URGnA+BGNFy6qnSzju78SuvSF9+6XtGQLSQEedDIEa0VKiQ2Axi4ECfswGih4w4HwIxoueii6Stt3bjhx6Sli71PSMgOri+lA+BGNGzzTYuGJu1a6Xhw33PCIgOri/lQyBGNNnytNWhNvfdJ23Y4HtGQDSQEedDIEY0NW0qtW/vxgsWSOPH+54REA1kxPkQiBFdW15l2rTJ52yAaCAjzodAjOiy4h777OPGU6ZIH3zge0ZA+JER50MgRnTRqxjIPK4v5UMgRrR17Chtv70bP/WU9MMPvmcEhBsFPfIhECPaqlSRevZ0Y2sCYSeoAaQPGXE+BGLg0ktdQDYjRki//cYzAdKFw1r5EIiBevWk885zz2HFCmnMGJ4JkC4c1sqHQAxseZVpwAB6FQPpQkacD4EYMC1bSkcd5Z7FnDnSiy/yXIB0Z8TlCUEEYqCwrPjee3k2QDoz4ooVeb6b8XYEiDnxRKlZMzeeNEn67DOeDZCujNhakiIPgRiIXybr2zdxrxhAapER51Nu0yYK7AJ/satLjRq509N2penHH5V3qhpAauy+uzRrlmtHav/OQEYMJKhRQ+ra1Y3Xr5eGDeMBAenIiFma/gtL08CWLrssOM05eLALyABSu0fMYa2/EIiBLe24o3TmmW68aJH0xBM8IyBVyIjzIRADyVxl4igFkBoc1sqHQAwU5KCDpP32c+NPPpHeeYfnBKQC15fyIRADBaFXMZAeZMT5EIiBwpx1ltSwoRs/95z07bc8K6CsyIjzIRADhalUyZ2gNrZHPGgQzwooKzLifCjoARRl6VKpcWNp7Vpp662l+fNdIQIApbPVVu7fkzVaoYxsHjJioCh16kidO7vxqlXSqFE8L6AsuL6UD4EYKE6fPsF44MDEfqoASoaCHvkQiIHi7LGHdNxxbvzdd9KECTwzoDTsrMXGjW5Micu/EIiBZNCrGCi7+NUkSlz+hUAMJMMyYusaY6y4x/TpPDegLIGYjPgvBGIg2QIf9CoGyoaMuEAEYiBZ558v1a7txo89Ji1cyLMDSnNQy5AR/4VADJTk/mP37m68YYM0dCjPDigJMuICEYiBkujZMzhkYoF43TqeH5AsMuICEYiBkmjUyNWgNkuWSI88wvMDkkVGXCACMVBS9CoGSoeMuEAEYqCk9t9fOvhgN/7iC+mNN3iGQDK4vlQgAjFQGhT4AMqWEVPQ4y8EYqA0Tj9datLEjZ9/Xpo9m+cIFIeMuEAEYqA07N18796JzSAAFI3DWgUiEAOldfHFUvXqbjx6tLRsGc8SKAqHtQpEIAZKa9ttpQsvdOM1a6SRI3mWQFHIiAtEIAZS1at40KDEd/wAEpERF4hADJTFLrtIJ5/sxvPmSc88w/MECkNGXCACMVBWXGUCkkNGXCACMVBWRx4ptWzpxpMnSx9/zDMFCsL1pQIRiIFU9ComKwaKR0GPAhGIgVQ491xpu+3ceNw4af58niuwJTLiAhGIgVSoWlXq0SN41z94MM8V2BIZcYEIxECqWCCuVMmNhw1zd4sBBMiIC0QgBlKlQQPpnHPc2KpsjR3LswXicX2pQARiIJX69g3G994rbdzI8wViuL5UIAIxkEqtW0uHHebGs2ZJr77K8wViyIgLRCAGUq1fv8SsGIBDRlwgAjGQaqecIu28sxu/8or05Zc8Y8CQEReIQAykWoUKic0gBgzgGQOGjLhABGIgHS66SNp6azd+6CFp6VKeM8D1pQIRiIF02GYb6eKL3XjdOmn4cJ4zQEGPAhGIgXTp3dvVoTb33Sf9/jvPGtFGRlwgAjGQLk2bSu3bu/GCBdL48TxrRBsZcYEIxECmrjL17y9t2sTzRnSREReIQAyk06GHSvvs48ZTp0offMDzRnRxfalABGIgk72KLSsGoorrSwUiEAPp1rGjtP32bvz009IPP/DMEU1kxAUiEAPpVqWK1KuXG1sTCDtBDUQRGXGBCMRAJnTv7gKyGTFC+u03njuih4y4QARiIBPq1ZPOO8+NV6yQHnyQ547oISMuEIEYyJT4Q1tWf5pexYgari8ViEAMZErLltJRR7nx3LnSiy/y7BEtFPQoEIEYyCSuMiHKyIgLRCAGMunEE6Xmzd34jTekzz7j+SM6OKxVIAIxkEnly0t9+wav6VWMKOGwVoEIxECmXXCBVLOmGz/yiLR4MV8DRAMZcYEIxECm1aghde3qxuvXS/ffz9cA0UBGXCACMeCrV3GFCm48ZIgLyEDYkREXiEAM+NCkiXTGGW68aJH0xBN8HRB+ZMQFIhAD2XKViV7FCDuuLxWIQAz4ctBB0n77ufGnn0rvvMPXAuFGQY8CEYgBn72K+/ULXt97L18LhBsZcYEIxIBPHTpIDRu68XPPSd98w9cD4UVGXCACMeBTpUrSZZe5se0RDxrE1wPhRUZcIAIx4Fu3blK1am48apS0cqXvGQHpwfWlAhGIAd/q1JE6d3bjVatcMAbCiOtLBSIQA9mgT59gPHBgYuYAhAUZcYEIxEA22GMP6bjj3Pi776QJE3zPCEg9MuICEYiBbMFVJkQpI7ZOZMjDkwCyxbHHSrvv7sZW3GP6dN8zAtKTEVuddbtHjzwEYiBb2Dem+F7FFPhAWDPiWMMT5Cm3aRMFboGssWaN1Lix9Ouv7o7xDz9IDRr4nhWQGi1bSl98IW21lbR6NU91MzJiIJvYN6ju3d14wwZp6FDfMwJSh4y4QARiINv07ClVrOjGFojXrvU9IyC1e8Sxv9/IQyAGsk2jRtJZZ7nxkiXSo4/6nhGQGmTEBSIQA9neq3jsWHoVI1yBmIw4AesDQDbaf39X9rJdO6lTJ9+zAVJ/fQl/IRAD2WrMGHdgy05PA2FARlwglqaBbEYQRpiQEReIQAyEyMaNG31PASgcGXGBCMRASLz99tvq37+/vv76a99TAQpGRlwgAjGQ42LF8VavXq3hw4fr5ptv9j0loGBcXyoQgRjIceU2F88/4YQTdPnll2vixImaTsMIZCMKehSIQAzkgIJKwtt+8Msvv6y5c+f+FZCPO+447b333rqXhhHIRmTEBSIQA1nOAm379u01f/78hKC8YcMGPf3003lZcOznGzdurD333FNffPGF1lgDCSBb2N9bDmsViEAMZLkPP/xQn3/+ed6P+EBcpUoVHXnkkXrxxRc1efLkvIy4UqVKmjVrlurVq/fXkjWQFeJP9FPQIwEFPYAst8cee2jVqlVavHhx3uvy5cvnLUvbf5s0aaLKlSurX79+atu2rdatW5cXlK+99lpVq1bN99SBQCwbNpS4TEAgBrJc69attf322+uVV17J2wOuX79+XhA2FnhtKbpBgwYaNmxYXoA+5phj1NnKYwLZeFDLkBEnYGkayAFXXXWVXnrpJY0YMSLv9a+//qqPP/5Y119/vTp06KBnn31Wr776at4y9gsvvKCmTZv6njKQiIy4UOU2FXQcE0DW6dSpkx5//PG8pWrLkOfMmaNly5ZpwoQJOvzww31PDyjasmVS7dpufPzx0ksv8cQ2IyMGcsQ999yT98MOZNmp6Hbt2uXtBxOEkXMZMUvTCciIgRyzdOlS1alTx/c0gJL5+WepQQM3bt9eevZZnuBmZMRAjiEIIyeREReKQAwAyOypaa4vJSAQAwDSj4y4UARiIAxlAzds8D0ToGhcXyoUgRjIZVZt67rrpJ12slNcvmcDFI6CHoUiEAO57M47pTvukBYskIYN8z0boHBkxIUiEAO5rFcv63/oxoMHS7//7ntGQMHIiAtFIAZymZWyPO00N7asePx43zMCCkZGXCgCMZDrNvcjztO/vzvABWQbMuJCEYiBXHfoodI++7jx1KnS5Mm+ZwTkx/WlQhGIgVxne8TxWfG99/qcDVAwCnoUikAMhEHHjlL9+m789NPSDz/4nhGQiIy4UARiIAyqVJF69nTjjRulQYN8zwhIREZcKAIxEBbdu7uAbEaOlFat8j0jIEBGXCgCMRAW9epJnTq58YoV0pgxvmcEBLi+VCgCMRAmffsG4wED3DI1kA24vlQoAjEQJi1bSkcd5cZz50ovvOB7RoBDRlwoAjEQNlxlQjYiIy4UgRgImxNPlJo3d+M33pA++8z3jAAy4iIQiIGwKV8+ca+YAh/IBmTEhSIQA2F0wQVSrVpu/Mgjrm8x4BPXlwpFIAbCqEYNqWtXN7bWiPff73tGiDoKehSKQAyE1WWXSRUquPGQIdL69b5nhCgjIy4UgRgIqyZNpDPOcONFi6THH/c9I0QZ15cKRSAGonSViV7F8IXDWoUiEANhdtBB0v77u/Gnn0rvvON7RogqMuJCEYiBKPUq7t/f52wQZWTEhSIQA2HXoYO0ww5uPGGC9M03vmeEKCIjLhSBGAi7SpXcCWpje8T0KoYPZMSFIhADUWB3iqtVc+NRo6SVK33PCFHD9aVCEYiBKKhTR+rc2Y1XrZIeeMD3jBA1FPQoFIEYiIr4+tMDByZmKEC6kREXikAMRMXuu0vHH+/G33/vDm4BmUJGXCgCMRAlXGWCL2TEhSIQA1Fy7LEuMzbvvitNm+Z7RogKri8VikAMRLnAx4ABPmeDKOH6UqEIxEDUdOok1a7txtYIYuFC3zNCFJARF4pADETNVltJ3bu78YYNrkUikG5kxIUiEANR1KuXVLGiG99/v7R2re8ZIezIiAtFIAaiyGpPn322Gy9ZIj36qO8ZIezIiAtFIAaiKr7Ah3Vlolcx0onrS4UiEANRZX2KDz7YjWfOlCZN8j0jhBkFPQpFIAairF+/YHzvvT5ngrAjIy4UgRiIstNOk5o0ceMXXpBmz/Y9I4QVGXGhCMRAlNnJ6d69g9cU+EC6kBEXikAMRN0ll0jVq7vxgw9Ky5b5nhHCiOtLhSIQA1FXq5bUpYsbr1kjjRzpe0YII64vFYpADEDq0yd4CoMGJX7TBFKBjLhQBGIAUvPm0sknuycxb5709NM8FaQWGXGhCMQAHK4yIZ3IiAtFIAbgHHGE1LKlG3/wgfTRRzwZpCcjLk/oicfTAFBwr2IKfCAdGbEFYfu7hr8QiAEEzj1XqlvXjceNk+bP5+kgtRlxrOsX/kIgBhCoWlXq0SPIYAYP5ukgtRlxhQo80S2U27SJlisA4vz8syt7uWGDtO227hR1rOAHUFp77il9+aW09dbSypU8xzhkxAAS1a8vnXOOG1uVrbFjeUJI3dI0GXE+BGIA+cUf2rL60xs38pSQmqVp9ojzIRADyG+ffaTDD3fjWbOkV1/lKaFsyIgLRSAGUHxW3L8/TwllQ0ZcKAIxgIKdcorUtKkbW0Y8cyZPCqVHRlwoAjGAgtmhmvhmEAMH8qRQelxfKhTXlwAUzq6ZNGokrVrl7hjbVabttuOJoeRq13an8K3ByOzZPME4ZMQACrfNNtLFF7vxunXS8OE8LZQOGXGhCMQAimbL07Ei/VZp6/ffeWIoOUpcFoqinwCKtvPOUvv20jPPSAsWSOPHu5rUQElcfbW0fn1Qyxx/YY8YQPHeeSe4V7zvvtLHH9NBByVjKynWdcl+UNQjAYEYQPGsJH2bNtInn7jX770ntW3LkwNSgD1iAMWzLKZfv+A1vYqBlCEjBpAc29/baSfXnckOb33zjXsNoEzIiAEkp0oVqWdPN7YmEPfdx5MDUoCMGEDyfvlFatzYZcc1a7oCH9ZfFkCpkREDSJ5dPenUyY1XrJDGjOHpAWVERgygZD7/XGrVyo2bNZO+/joo+AFs4c/NFbXKlSun8vw9KRCBGEDJHX20NGmSG0+Y4Do1AZstW7ZM33zzjT7//HMtWLBA69evV7Vq1bTDDjuodevW2mWXXVS5cmWe12YEYgAl98IL0sknu/GRRwZBGZE3ceJEXXfddfriiy/ygm3FzcU7LDO2H3/88UdeIL7zzjt1yimn5GXKUcd6EoCSO+EE10XHvPGGNGMGTxF5Qfi0005T1apVNWzYMH344YeaO3eu5s2bl5chf/DBBxo6dKi22WabvN833sqlgowYQClZA4jLLnPjLl2kUaN4lBHXpk0bbbvttnr99deL/b0nn3yyvv/+e02dOjUvcEcZGTGA0rngAqlWLTd+5BFp8WKeZMTNmTNHZ5xxxl+vN1lp1EKcd955mj17tjbanfSIIxADKJ0aNaSuXYOC/vffz5OMuCZNmui9997T6tWr814Xtf9re8iWPVekAQSBGEAZ2NJ0hQpuPGSIK/SByOrVq5fGjRunq666SjNnztTSpUu1atUqrVmzRmvXrs37sXjxYt1xxx266667dP7556tSpUqKOk5NAyibs8+Wxo1z4wcfdEvWiCw7MW0HspYvX66GDRtq++23z9sDttPSS5Ys0XfffZd3n7hjx44aNWqUqljp1IgjEAMomw8+kA4+2I333luaPp1exRG3cOFCvfrqq/r000/1008/aeXKlXlL0HXq1NHuu++ed1CrRYsWvqeZNQjEAMrGDuQceKD08cfu9ZtvSu3a8VSBJHFYC0DZ2IGcyy8PXtOrGCgRMmIAZbdhg7TzztJPP7nAPGeO9Le/8WSBJJARAyg7O/kaK+5hS9WDBvFUkcj+XtidYe4N50NGDCA1li51vYrXrnV3jOfPdz2LgS+/lAYMcFfdTjwxqFOOPGTEAFKjTh2pc2c3/u03Sl4iMG+eNHy4NHSoNHUqT2YLBGIAqdO3bzAeONBa7vB0If3xR/AUYgVg8BcCMYDU2X136fjj3fj776XnnuPpIvENGSUt8yEQA0gtrjJhS2TERSIQA0itY491mbF5911p2jSecNTFZ8QsTedDIAaQ3gIfdloW0RafEbM0nQ+BGEDqdeok1a7txo8/bsWHecpRRkZcJAIxgNTbaiupe/eg6pa1SER0kREXiUAMID169QqWIe+/3xX6QDSREReJQAwgPXbYwfUqNkuWSI88wpOOKq4vFYlADCBzV5ms3jCih+tLRSIQA0if/faT2rZ145kzpUmTeNpRREZcJAIxgPSiwAfIiItEIAaQXqedJjVp4sYvvCB9/TVPPGrIiItEIAaQXnZyunfvxGYQiBYy4iIRiAGk3yWXSNWru/GDD0rLlvHUo4TrS0UiEANIv1q1pC5d3HjNGmnECJ56lFDQo0gEYgCZ0aePq0NtBg1yFbcQDWTERSIQA8iM5s2lk0924/nzpWee4clHBRlxkQjEADKHq0zRREZcJAIxgMw54gipZUs3/uAD6aOPePpRwPWlIhGIAfjrVWxlLxF+XF8qEoEYQGade65Ut64bjxvn9osRbmTERSIQA8isqlWlHj2Cb9D33cdXIOzIiItEIAaQeRaIK1d24+HDpdWr+SqEGRlxkQjEADKvfn3pnHPc2KpsjR3LVyHMyIiLRCAG4EffvsF4wABp40a+EmHF9aUiEYgB+LHPPtLhh7vxrFnSK6/wlQgrCnoUiUAMwB+uMkUDGXGRCMQA/DnlFKlpUzd+9VVp5ky+GmHEYa0iEYgB+FOhgmsGEb9XjPDhsFaRCMQA/LL2iFtv7cZ2enrJEr4iYUNGXCQCMQC/ttlGuuQSN163zt0rRriQEReJQAzAv969pfKbvx0NHiz9/rvvGSGVyIiLRCAG4N/OO0vt27vxggWuBjXCg4y4SARiANl5lWnTJp+zQSpxfalIBGIA2eHQQ6XWrd146lRp8mTfM0KqUNCjSARiANnZq7h/f5+zQSqREReJQAwge3Ts6BpCmGeekb7/3veMkMqM2A7k2RsuJCAQA8ge1hqxVy83tiYQ9CoOV0ZsBVyQT7lNmzgRASCL/PKL1LixtH69VLOmNG9eUPADuWmvvaTPPpOqVZPWrPE9m6xDRgwgu9StK3Xq5MYrVkgPPuh7RkjV0jQZcYEIxACyD72Kw7k0XbGi75lkJQIxgOzTsqV09NFu/M030gsv+J4RyoKMuEgEYgDZiatM4UFGXCQCMYDsdMIJ0i67uPGbb0ozZvieEUqLjLhIBGIA2cnunG65V4zcxPWlInF9CUD2+u03d5Vp+XJ3x9iuMtWr53tWKCn7mtm1NGvu8e23PL8tkBEDyF41akhdu7qxtUYcOtT3jFAaZMRFIhADyG6XXRbcPx0yxBX6QG7uEXN9qUAEYgDZrUkT6cwz3XjxYunxx33PCCVFRlwkAjGA3LvKRGXe3ML1pSIRiAFkvwMPlPbf343tGtPbb/ueEUqC60tFIhADyH7WOq9fv+D1vff6nA1Kioy4SARiALnB9ol32MGNJ0xwpS+R/aydZWwrgaYPBSIQA8gNlSq5E9TGvrEPHOh7RihJNmw4NV0gAjGA3NGtm+tpa0aNcm0SkRv7w4aMuEAEYgC5o3Zt6YILgqpbFoyROxkxgbhABGIAuaVPn2Bsy9Px3+iR3RkxS9MFIhADyC277y4df7wbf/+99NxzvmeEopARF4tADCD3cJUpd3BYq1gEYgC555hjXGZs3n1XmjbN94xQGA5rFYtADCA3C3zEl72kwEf2IiMuFoEYQG7q1MmdojZPPCEtWOB7RigIGXGxCMQActNWW0mXXurGGzbQqzhbkREXi0AMIHf17Blcibn/fmntWt8zwpbIiItFIAaQu6z29Nlnu/GSJdIjj/ieEbbE9aViEYgBhOsqE72KswsFPYpFIAaQ2/bdV2rb1o1nzpQmTfI9I8QjIy4WgRhA7ou/ytS/v8+ZYEtkxMUiEAPIfaedJu24oxu/+KL09de+Z4QYMuJiEYgB5D47Od27d/CaXsXZg+tLxSIQAwiHiy+Wqld34wcflH791feMYLi+VCwCMYBwqFVL6tLFjdeskUaO9D0jGDLiYhGIAYSrV7HVoTaDBrmKW/CLjLhYBGIA4dG8uXTyyW48f770zDO+ZwQy4mIRiAGEC1eZsgsZcbEIxADC5YgjpFat3PjDD90P+MP1pWIRiAGEu1fxgAE+ZwMKehSLQAwgfM45R6pb143HjZPmzfM9o+giIy4WgRhA+FStKvXoEQSCwYN9zyi6yIiLRSAGEE4WiCtXduPhw6XVq33PKJrIiItFIAYQTvXruyVqs2yZNHas7xlFE9eXikUgBhBeffsm9ireuNHnbKKJ60vFIhADCK999pEOP9yNrSPTK6/4nlH0kBEXi0AMINz69UvMipFZZMTFIhADCDcredm0qRu/+qo0c6bvGUULGXGxCMQAwq1CBdcMIoYCH5lFRlwsAjGA8LP2iFtv7cZ2enrJEt8zig6uLxWLQAwg/LbZRrrkEjdet87dK0ZmUNCjWARiANHQu7dUfvO3vPvuk37/3feMooGMuFgEYgDRsPPOUvv2brxwoatBjfTjsFaxCMQAonuVadMmn7OJBg5rFYtADCA6DjlEat3ajadOld5/3/eMwo+MuFgEYgDR7VVMgY/0IyMuFoEYQLR07OgaQphnnpG+/973jMKNjLhYBGIA0WKtEXv1cmNrAmEnqJE+ZMTFIhADiJ7u3aUqVdx4xAhp1SrfMwovri8Vi0AMIHrq1pU6dXLjlSulBx/0PaPwoqBHsQjEAKIp/tCW1Z+mV3F6kBEXi0AMIJpatJCOPtqNv/lGev553zMKJzLiYhGIAUQXV5nSj4y4WARiANF1wgnSLru48ZtvSjNm+J5R+HB9qVgEYgDRZU0g+vYNXtOrOPW4vlQsAjGAaOvcWapVy40feURatMj3jMKFjLhYBGIA0VajhtS1qxtba8T77/c9o3AhIy4WgRgALrtMqlDBPYchQ6T163kmqUJGXCwCMQA0aSKdeaZ7DosXS489xjNJFTLiYhGIAaCgq0z0Kk4Nri8Vi0AMAObAA6X993fPwq4xvf02zyUVKOhRLAIxAMR6FffrFzwLehWnBhlxsQjEABBj+8Q77ODGEyZIc+fybMqKPeJiEYgBIKZSJeWdoDa2RzxoEM8mVRmxrThYARXkU27TJk4kAMBffv1VatRIWrvW3TGeP1+qWZMHVFpt2kjTp7s3OXZPG/nw9gQA4tWuLV1wgRv/9pv0wAM8n1QsTcfuaSMfAjEAbKlPn2A8cGDiPidKtzRdsSJPrhAEYgDY0u67S8cf78Y//OAObqF0yIiLRSAGgIJwlSk1yIiLRSAGgIIcc4zLjM2770rTpvGcSoOMuFgEYgAoiF232bLsJUqfEXNYq1BcXwKAwtgVpsaNpaVL3fWb77+XGjbkeZWEPa+FC91z/PFHnl0ByIgBoDDVqkndu7vxhg2uRSJKhoy4WARiAChKz57B1Zv773dZMpLHYa1iEYgBoChWe/rss93YlqgfeYTnVRIc1ioWgRgASnqVicrAySMjLhaBGACKs+++Utu2bjxzpvT66zyzZJERF4tADADJ4CpT6ZARF4tADADJOO00accd3fjFF6Wvv+a5JYOMuFgEYgBIhp2c7t07eD1gAM+tOBs3BvvpFPQoFIEYAJJ18cVS9epuPGaM612M4pelDd2XCkUgBoBk1aoldenixmvWSCNH8uySDcRkxIUiEANASXsVWx1qM2iQq7iFgsX3cSYjLhSBGABKonlz6eST3Xj+fOnpp3l+hSEjTgqBGABKiqtMyWGPOCkEYgAoqSOOkFq1cuMPP3Q/UPTSNHvEhSIQA0BJ0as4OWTESSEQA0BpnHOOVLeuG48fL82bx3PcEhlxUgjEAFAaVatKPXoEmd/gwTzHLZERJ4VADAClZYG4cmU3Hj5cWr2aZxmPjDgpBGIAKK369d0StVm2THroIZ5lPK4vJYVADACpuspk9aetvjIcCnokhUAMAGWx995Su3ZubB2ZXnmF5xlDRpwUAjEApDIr7t+f5xlDRpwUAjEAlJWVvGza1I1fe02aOZNnasiIk0IgBoCysqpR1gwihl7FDteXkkIgBoBUuOgiaZtt3HjsWGnJEp4r15eSQiAGgFTYemvp4ovdeN06adgwnisZcVIIxACQKr17S+U3f1u1Slu//x7tZ0tGnBQCMQCkys47S6ed5sYLF0rjxkX72ZIRJ4VADADpvMq0aVN0ny8ZcVIIxACQSoccIrVu7cbTpknvvx/d58v1paQQiAEglehVHKCgR1IIxACQah07uoYQ5plnpO+/j+YzJiNOCoEYAFLNWiP26uXG1gRi0KBoPmMOayWFQAwA6dC9u1SlihuPHCmtWhW958xhraQQiAEgHerWlc4/341XrpQefDB6z5mMOCkEYgBIl759E+tPxwemKCAjTgqBGADSpUUL6eij3fibb6QXXojWsyYjTgqBGAAyVeDj3nuj9azJiJNCIAaAdDrhBGmXXdz4zTelGTOi87y5vpQUAjEApJM1gYjfK45SVkxBj6QQiAEg3Tp3lmrVcuNHH5UWLYrGMycjTgqBGADSrUYNqVs3N7bWiPffH41nTkacFAIxAGSCVdqqUMGNhwyR1q0L/3MnI04KgRgAMqFJE+nMM9148WLp8cfD/9y5vpQUAjEA+LrKFPZexVxfSgqBGAAy5aCDpAMOcGO7xvT22+F+9mTESSEQA4CvrLh//3A/ezLipBCIASCTbJ+4USM3njhRmjs3vM+fjDgpBGIAyKRKlaTLLnNj2yMOc69iMuKkEIgBINO6dpWqVXPjUaOkFSvC+TXg+lJSCMQAkGm1a0sXXODGv/0mPfBAOL8GFPRICoEYAHyIrz89cGBi0AoLMuKkEIgBwIfddnOdmcwPP0gTJoTv60BGnBQCMQD4EvarTGTESSEQA4Avxxwj7bGHG7/3njR1ari+FlxfSgqBGAB8KVcuMSseMCBcXwuuLyWFQAwAPnXqJNWp48bWCGLBgvB8PciIk0IgBgCf7D5x9+5BBmktEsOCjDgpBGIA8K1nT6liRTe+/35p7VqFAhlxUgjEAODbDjtIHTu68dKl0iOPKBTIiJNCIAaAbBDGXsVcX0oKgRgAssG++0pt27rxzJnS668r51HQIykEYgDIFv36JWbFuY6MOCkEYgDIFu3bSzvu6MYvvijNmqWcxmGtpBCIASBb2Mnp3r0Tm0HkMg5rJaXcpk1hOBEAACGxfLnUqJG0erW01VbSvHmubWIuOuAA6eOPXQWxjRt9zyZrkREDQDapVUu66CI3XrNGGjlSOZ8RV6jgeyZZjUAMANnGlqctizSDBkkbNiin94hjxUpQIAIxAGSb5s2lk0924/nzpaefVk4iI04KgRgAslEYrjLFMmKWpotEIAaAbNSundSqlRt/+KH7kasZMUvTRSIQA0Au9CrOxayYjDgpBGIAyFbnnCPVq+fG48e7q0y5hIw4KQRiAMhWVatKPXoE2eXgwcopZMRJIRADQDa79FKpcmU3Hj7cFfrIFVxfSgqBGACyWf36bonaLFsmPfSQcgbXl5JCIAaAbBd/aGvAgNwpF0lGnBQCMQBku733dteZzNdfSy+/rJxARpwUAjEA5IJcvMpERpwUAjEA5AIredm0qRu/9po0c6ayHhlxUgjEAJALrExk376Je8XZjutLSaEfMQDkilWrXK/ilSvdHWMr8LHddspK1uq+/OZc78ADpQ8+8D2jrEVGDAC5YuutpYsvduN166Rhw5T12bCh6UORCMQAkGu9imOZplXa+v13ZfX+sKHpQ5EIxACQS3beWTrtNDdeuFB68kllJTLipBGIASCXrzL17+/2Y7M5EJMRF4lADAC55pBDpNat3Xj6dOn995XVS9PsEReJQAwAudiruF+/7C7wQUacNAIxAOSis892DSHMM89I332nrEJGnDQCMQDkImuN2KuXG1sTiPvuU1YhI04agRgAclX37q6whxk50hX8yBZkxEkjEANArqpbV+rUyY2t2tbo0coaXF9KGoEYAHJZfP3pgQMTA6BPFPRIGoEYAHJZixbS0Ue78TffSC+8oKxARpw0AjEA5LpsvMrEYa2kEYgBINcdf7y0yy5u/Oab0qef+p4Rh7VKgEAMALnOmkBkW69iMuKkEYgBIAwuuECqVcuNH31UWrTI73y4vpQ0AjEAhEH16lK3bm5srRGHDvU7HzLipBGIASAsrNJWrMGCBeJ16/zNhYw4aQRiAAiLJk2kM89048WLpccf9zcXri8ljUAMAGG+yuSrVzEFPZJGIAaAMDnwQOmAA9x4xgzprbf8zIOMOGkEYgAIm8sv91/gg4w4aQRiAAgb2ydu1MiNJ06U5s7N/BzIiJNGIAaAsKlUSbrsMje2PWJrBpFpXF9KGoEYAMKoa1epWjU3HjVKWr48s5+f60tJIxADQBjVru2qbZnVq10wziQy4qQRiAEgKr2K47PUdCMjThqBGADCarfdpBNOcOMffpCeey5zn5uMOGkEYgAIM19XmciIk0YgBoAwO+YYaY893Pi996SpUzPzebm+lDQCMQCEWblyfrJiCnokjUAMAGHXqZNUp44bP/GEtGBB+j8nGXHSCMQAEHZ2n7h79yBTHTIk/Z+TjDhpBGIAiEqv4ooV3fj++6W1a9P7+ciIk0YgBoAoaNhQ6tjRjZculR5+OL2fj+tLSSMQA0BUbHloK529irm+lDQCMQBExb77Socc4sZffim9/nr6PhcZcdIIxAAQ1ay4f//0fR4y4qQRiAEgStq3l3bc0Y1fekmaNSs9n4eMOGkEYgCIEjs53bt38DpdvYrJiJNGIAaAqLnkEqlGDTceM0b69dfUfw6uLyWNQAwAUVOzptSlixuvWSONGJH6z0FBj6QRiAEgivr0cXWozX33SRs2pPbjkxEnjUAMAFHUrJl0yiluPH++9PTTqf34HNZKGoEYAKIqnVeZOKyVNAIxAERVu3ZSq1Zu/NFH0ocfpu5jkxEnjUAMAFGVzl7FZMRJIxADQJSdc45Ur54bjx8vzZuXmo9LRpw0AjEARFnVqlKPHkHwtBPUqUBGnDQCMQBEnQXiypXdePhwafXqsn9Mri8ljUAMAFG3/fbSuee68fLl0kMPlf1jUtAjaQRiAIDUt2/ioa2NG8v2VMiIk0YgBgBIe+/trjOZ2bOll18u21MhI04agRgA4KTyKhMZcdIIxAAA5+STpb/9zY1fe02aObP0T4brS0kjEAMAnAoVXDOIVGTFXF9KGoEYABCw9ojbbOPGY8dKv/xSuqdDRpw0AjEAILD11tIll7jx+vXuXnFZM+LyhJqi8HQAAIkuuywInoMHS7//XvqMuGJFnm4xCMQAgEQ77yyddpobL1woPflk6TNi23dGkQjEAIDiexVv2lS6jJhAXCwCMQAgv0MOkdq0cePp06X33y9dRszSdLEIxACA4nsVW1ZcEmTESSMQAwAKdvbZUoMGbvzss9J33yX/pMiIk0YgBgAUzFoj9urlxtYEoiS9ismIk1Zu06aS7sADACLDCno0aSKtW+cKfcyf7+4aF6dRI+mnn9x/583LxExzFhkxAKBwdetKnTq58cqV0ujRyT0tri8ljUAMAEi+V/HAgYnlKwtDQY+kEYgBAEVr0UI65hg3/uYb6YUXin9iZMRJIxADAIpX0qtMZMRJIxADAIp3/PHSrru68VtvSZ9+WvTvJyNOGoEYAJBEtCifuFc8YEDRv5/rS0nj+hIAIDmrV7vrSMuXuzvGP/wg1a9f8O+1GtN293jffaUpU3jCRSAjBgAkp3p1qVs3N7bWiPffX/Dvs/IUFoQNTR+KRSAGAJSsV3EsuA4d6gp9bCkWhA1NH4pFIAYAJK9xY6lDBzdevFh6/PHCD2oZMuJiEYgBAKntVRxf8IOMuFgVi/8tAADEOfBA6YADpI8+kj77zF1nOuKIvHi8YIE0/T1psS7S76qsyj83V72JUuvWUsOGrrsiEnFqGgBQck88If3979qocnr9wBs0otHNevddadGiwv/I9ttLhx4qde0qHX20uxEFAjEAoBT+XLdBwxvcqHuWd9FcNS/xn2/WTLriCncIO+rbyGTEAIASmTNH6tJFev/9xJ+vVctdG27TRvrb39xVY7vlZOWpp02Tpk51V5DjtW3rGjo1L3ksDw0CMQAgaWPHuiw2/tbSUUdJPXtKp55a9NksO0w9YYI0ZIg0aVLw81WrSsOHS+efH80vBIEYAJAUC6C9egWvLesdNUo67LCSP8B33nFZ9bffBj83eLAL6FHDVjkAIKlMOD4I24GrGTNKF4SN/Tk7cN21a/Bz9vHt80QNGTEAoNg94VatguXoa66Rbr01NVeRNm2SrrtOuu0297paNRfgo7RnTCAGABTKanMcfnhwMMsy2GHDUnsfeNMmqXt3acSI4ADX229H5zQ1gRgAUCgrJx3bt7U9YctWrfdDOho7tWoV7BnbfnSPHtH4wrBHDAAokPVuuOee4LUdzCouCP/222+68cYbdfzxx6t27doqV66cHnzwwWKfcPXq7uPHV86M7x0RZgRiAECBXn9dmjs3uKKUzMGsJUuW6N///re++uor7bXXXiV6socfLh15ZLAvHX/FKcwIxACAAsX2bE2y14oaNGighQsX6ocfftBdd91V4ifbM+7z2N3iKCAQAwAKPEBltaNjFbOsWEcyqlSpovr165f6ibZv7z6fee+9aHxhCMQAgHysi1KsgYOVrcxUN8OKFV2JTPPzz24eYUcgBgDkM316MI4FxkxpE/f5rEZ12BGIAQD5LF4cjO3aUib9Le7z/fKLQo9ADADIx7omxVgXpUyqUiUYr1+v0CMQAwDyiQ++8UE5E9avLzgohxWBGACQT716wdj6CWfSN3Gfr25dhR6BGACQT+vW/g5MTZvm76CYDxk6kA4AyCUNG0rbb++uME2dKv3xR/JXmO677z4tX75cCzbfPZo4caLmz5+fN+7du7dq1qxZ6J/9448gENt1ZJtH2NH0AQBQoLPOksaPd+OnnpLOOCO5B7XTTjvlVdYqyHfffZf364V56impQwc3tv+OGxf+Lw5L0wCAAlnLwxjrhpSs77//Xps2bSrwR1FBeMvP062bIoGMGABQIOt+tOuuQeMH6xGcTOOH0nr7baldOzdu3lyaNUsqH4F0MQL/iwCA0rAgeMUVweuLLnJ9g9Nh9Wr38WP69YtGEDYR+d8EAJSGLQ+3bRtcK7IAaQ0hUmnTJunyy6Vvv3WvDzlE6t5dkUEgBgAUqkIFafRoqWrVoDXiddelLhhv2iRde600cqR7Xa2aNGpUdLJhE6H/VQBAadh+bXxv4NtucxlrWZep7c936ybdfnvwc8OGuc8XJRzWAgAkxU409+oVvG7a1GWvhx9euoNZtif87eblaDN4sNSzZ/S+GARiAEDSxo512fDatcHPHXmkC6Dt2xdd9GPDBmnCBBfQ33gj+Plq1VwmfP750fxCEIgBACUyZ47UpYv0/vuJP1+rlitJaT+slaE1bLAGDnbIy6pl2Y/lyxP/TNu2bg86asvR8QjEAIAS+/NPt298zz3BPeOSaNbMXY3q1s0dCIsyAjEAoExFPyZNckH5vfekn38u/Pda7Wi7mmTB96ijonUyuigEYgBAylifB1uC/uUXtyxty9PWytCWq6PQwKE0CMQAAHjEwgAAAB4RiAEA8IhADACARwRiAAA8IhADAOARgRgAAI8IxAAAeEQgBgDAIwIxAAAeEYgBAPCIQAwAgEcEYgAAPCIQAwDgEYEYAACPCMQAAHhEIAYAwCMCMQAAHhGIAQDwiEAMAIBHBGIAADwiEAMA4BGBGAAAjwjEAAB4RCAGAMAjAjEAAB4RiAEA8IhADACARwRiAAA8IhADAOARgRgAAI8IxAAAeEQgBgDAIwIxAAAeEYgBAPCIQAwAgEcEYgAAPCIQAwDgEYEYAACPCMQAAHhEIAYAwCMCMQAAHhGIAQDwiEAMAIBHBGIAADwiEAMA4BGBGAAAjwjEAAB4RCAGAMAjAjEAAB4RiAEA8IhADACARwRiAAA8IhADAOARgRgAAI8IxAAAeEQgBgDAIwIxAAAeEYgBAPCIQAwAgEcEYgAAPCIQAwDgEYEYAACPCMQAAHhEIAYAwCMCMQAAHhGIAQDwiEAMAIBHBGIAADwiEAMA4BGBGAAAjwjEAAB4RCAGAMAjAjEAAB4RiAEA8IhADACARwRiAAA8IhADAOARgRgAAI8IxAAAeEQgBgDAIwIxAAAeEYgBAPCIQAwAgEcEYgAAPCIQAwDgEYEYAACPCMQAAHhEIAYAwCMCMQAAHhGIAQDwiEAMAIBHBGIAADwiEAMA4BGBGAAAjwjEAAB4RCAGAMAjAjEAAB4RiAEA8IhADACARwRiAAA8IhADAOARgRgAAI8IxAAAeEQgBgDAIwIxAAAeEYgBAPCIQAwAgEcEYgAAPCIQAwDgEYEYAACPCMQAAHhEIAYAwCMCMQAAHhGIAQDwiEAMAIBHBGIAADwiEAMA4BGBGAAAjwjEAAB4RCAGAMAjAjEAAB4RiAEA8IhADACARwRiAAA8IhADAOARgRgAAI8IxAAAeEQgBgDAIwIxAAAeEYgBAPCIQAwAgEcEYgAAPCIQAwDgEYEYAACPCMQAAHhEIAYAwCMCMQAAHhGIAQDwiEAMAIBHBGIAADwiEAMA4BGBGAAAjwjEAAB4RCAGAMAjAjEAAB4RiAEA8IhADACARwRiAAA8IhADAOARgRgAAPnz/0kDVEdR+YuMAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 500x1415.35 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "np.random.seed(0)\n",
    "dn = op.network.Delaunay(points=4, shape=[1, 1, 0])\n",
    "op.visualization.plot_tutorial(dn);"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The basic premise of how OpenPNM stores topology can be stated in 1 sentence:\n",
    "\n",
    "**The pores on either end of a throat are just another property to be stored, along with diameter, length, etc.**\n",
    "\n",
    "In other words, referring to the above diagram, throat No. 5 has pores 1 and 4 on its ends.  Using the spreadsheet analogy, this implies a new column that stores the pair of pores connected by each throat.  OpenPNM calls this property ``'throat.conns'``: "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2021-06-24T11:25:06.713161Z",
     "iopub.status.busy": "2021-06-24T11:25:06.712287Z",
     "iopub.status.idle": "2021-06-24T11:25:06.715218Z",
     "shell.execute_reply": "2021-06-24T11:25:06.715912Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[0 1]\n",
      " [0 2]\n",
      " [0 3]\n",
      " [1 3]\n",
      " [2 3]]\n"
     ]
    }
   ],
   "source": [
    "print(dn['throat.conns'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Inspection of this *Nt-by-2* array shows that each row contains the pore indices corresponding to that throat."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "```{note} \n",
    "  **Topology Rules**: OpenPNM requires that the following be true about the network:\n",
    "  1. A throat is a connection between exactly two pores, no more and no less\n",
    "  2. Throats are non-directional, meaning that flow in either direction is equal\n",
    "  Other general, but non-essential rules are:\n",
    "  3. Pores can have an arbitrary number of throats, including zero; however, pores with zero throats lead to singular matrices and other problems so should be avoided.\n",
    "  4. Two pores are generally connected by no more than one throat. It is technically possible in OpenPNM to have multiple throats between a pair of pores, but it is not rigorously  supported so unintended results may arise. For instance, when converting to CSR format the dual throats are merged into one, so re-converting back results in loss of such throats.\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Adjacency Matrices: Dense vs Sparse\n",
    "\n",
    "The topology storage scheme described above is actually an [adjacency matrix](http://en.wikipedia.org/wiki/Adjacency_matrix), in a [sparse storage format](http://en.wikipedia.org/wiki/Sparse_matrix) known as IJV or COO.  \n",
    "\n",
    "An adjacency matrix is a *Np*-by-*Np* matrix with non-zero values at location (*i*, *j*) indicating that pores *i* and *j* are connected.  "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "OpenPNM Network objects have a method that generates the sparse adjacency matrix:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "<COOrdinate sparse matrix of dtype 'int64'\n",
      "\twith 5 stored elements and shape (4, 4)>\n",
      "  Coords\tValues\n",
      "  (0, 1)\t0\n",
      "  (0, 2)\t1\n",
      "  (0, 3)\t2\n",
      "  (1, 3)\t3\n",
      "  (2, 3)\t4\n"
     ]
    }
   ],
   "source": [
    "am = dn.create_adjacency_matrix(weights=dn.Ts, triu=True)\n",
    "print(am)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This can be converted to a *dense* array and visualized:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[0 1 1 1]\n",
      " [1 0 0 1]\n",
      " [1 0 0 1]\n",
      " [1 1 1 0]]\n"
     ]
    }
   ],
   "source": [
    "am = dn.create_adjacency_matrix()\n",
    "print(am.todense())"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In dense format the *Np-by-Np* shape is evident. The above matrix is symmetrical, so connections between pore *i* and *j* as well as between *j* and *i* are shown."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "An important feature of the adjacency matrix is that it is highly sparse (mostly zeros) and can be stored with a variety of [sparse storage schemes](https://docs.scipy.org/doc/scipy/reference/sparse.html) offered by Scipy.  OpenPNM stores the adjacency matrix in the 'COO' format, which essentially stores the coordinates (I,J) of the nonzero elements in an two-column wide array.\n",
    "\n",
    "Thus the throat property called ``'throat.conns'`` is an *Nt*-by-2 array that gives the index of the two pores on either end of a given throat.  "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "```{note} \n",
    "\n",
    "  **Additional Thoughts on Sparse Storage:**\n",
    "  * In pore networks there is (usually) no difference between traversing from pore *i* to pore *j* or from pore *j* to pore *i*, so a 1 is also found at location (*j*, *i*) and the matrix is symmetrical.\n",
    "  * Since the adjacency matrix is symmetric, it is redundant to store the entire matrix when only the upper triangular part is necessary.  The ``'throat.conns'`` array only stores the upper triangular information, and *i* is always less than *j*.\n",
    "  * Although this storage scheme is widely known as *IJV*, the ``scipy.sparse`` module calls this the Coordinate or *COO* storage scheme.\n",
    "  * Some tasks are best performed on other types of storages scheme, such as *CSR* or *LIL*.  OpenPNM converts between these internally as necessary, but users can generate a desired format using the ``create_adjacency_matrix`` method which accepts the storage type as an argument (i.e. ``'csr'``, ``'lil'``, etc).  For a discussion of sparse storage schemes and the respective merits, see this [Wikipedia article](http://en.wikipedia.org/wiki/Sparse_matrix).\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## ``conduit`` data\n",
    "\n",
    "A *conduit* in OpenPNM refers to a 'pore-throat-pore' connection, where the pores and throats are each called *elements*. These *elements* act as resistors (or conductors) in series. Each *conduit* therefor contains exactly 1 throat and 2 pores. Transport within the network occurs via these interconnected *conduits*, so the geometric properties of these elements are often needed for computing things like diffusive conductance. This is common enough that OpenPNM V3 includes a new helper function for obtaining this data, as illustrated below:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[0.47440675 0.23720338 0.55759468]\n",
      " [0.55759468 0.25069084 0.50138169]\n",
      " [0.50138169 0.2362208  0.47244159]]\n"
     ]
    }
   ],
   "source": [
    "D = pn.get_conduit_data('diameter')\n",
    "print(D)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The above array contains 3 columns, with the center column representing the throat diameter:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[0.23720338 0.25069084 0.2362208 ]\n"
     ]
    }
   ],
   "source": [
    "print(pn['throat.diameter'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "And the left and right columns represent the diameter of the pores on each end of the throat. The order of \"left\" vs \"right\" reflects the pore indices for the ``'throat.conns'`` array.  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[0 1]\n",
      " [1 2]\n",
      " [2 3]]\n"
     ]
    }
   ],
   "source": [
    "print(pn.conns)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In fact, the values for the pore elements can be obtained using numpy's fancy indexing:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[[0.47440675 0.55759468]\n",
      " [0.55759468 0.50138169]\n",
      " [0.50138169 0.47244159]]\n"
     ]
    }
   ],
   "source": [
    "R1_R2 = pn['pore.diameter'][pn.conns]\n",
    "print(R1_R2)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Which can be seen to match the diameter values in ``D`` above. "
   ]
  }
 ],
 "metadata": {
  "@webio": {
   "lastCommId": null,
   "lastKernelId": null
  },
  "kernelspec": {
   "display_name": "openpnm",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.13.0"
  },
  "toc": {
   "base_numbering": 1,
   "nav_menu": {},
   "number_sections": false,
   "sideBar": true,
   "skip_h1_title": true,
   "title_cell": "",
   "title_sidebar": "Contents",
   "toc_cell": true,
   "toc_position": {},
   "toc_section_display": true,
   "toc_window_display": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}
